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INTRODUCTORY  SUMMARY 


This  report  considers  the  errors  attributable  to  the  effects  of  non-ideal 

filter  transmission  characteristics  (e.g.,  finite  slopes  for  the  filter  skirts)  and 

varying  slopes  of  the  spectrum  applied  to  the  filter  input,  upon  the  measurement, 

correction,  or  extrapolation  (in  the  presence  of  masking  by  background  noise)  of 

aircraft  noise  data.  This  report  is  the  third  in  a  series  of  reports  on  correction 
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procedures  for  the  evaluation  of  aircraft  noise. 

Current  engineering  practice  for  evaluation  of  the  spectral  content  of 
aircraft  noise  is  well  defined  in  FAA  Part  36,^  and  in  a  number  of  related  standards 
and  documents  covering  the  details  of  aircraft  noise  measurement  and  spectrum 
analysis.^"®  In  current  engineering  practice,  exemplified  by  the  procedures  cited  in 
these  references,  no  corrections  are  applied  for  filter  sideband  or  spectrum  slope 
effects  when  analyzing  the  aircraft  spectrum  data. 

The  basic  approach  utilized  in  this  report  for  analyzing  these  errors  due  to 
the  filter  skirts  and  signal  spectrum  slopes  involves,  first,  computing  a  close 
approximation  to  the  true  spectrum  shape  of  the  signal  at  all  frequencies.  This  is 
one  of  the  key  departures  from  recent  related  studies  which  have  focused  on  just 
one  aspect  of  the  problem  -  specifying  attenuation  of  a  band  of  noise  due  to 
atmospheric  absorption,  given  a  definition  of  attenuation  at  single  frequencies.^  ^ 
The  method  outlined  in  this  report  is  applicable  to  this  and  other  problems 
encountered  in  evaluation  or  correction  of  aircraft  spectra.  The  method  is  based 
on  the  use  of  a  curve-fitting  spline  function  to  estimate  the  spectrum  shape 
between  center  frequencies  of  three  adjacent  one-third  octave  bands.  The  method 
includes  a  critical  iteration  routine  to  ensure  that  the  integration  of  this  estimated 
spectrum  over  the  effective  transmission  bandwidth  of  each  one-third  octave  band 
filter  produces  the  same  band  level  as  actually  specified  by  the  measured  value.  It 
is  impossible,  of  course,  to  know  the  exact  shape  of  the  aircraft  noise  spectrum  at 
all  frequencies  between  the  nominal  band  center  frequencies,  since  only  the 
integrated  band  level  is  specified  by  the  measured  data.  However,  the  analytical 
methods  employed  here  for  estimating  and  verifying  this  spectrum  shape  provide 
one  rational  and  consistent  approach  to  the  problem.  Although  the  precise 
quantitative  results  obtained  with  this  iterative  technique  are  dependent  upon  the 
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exact  form  of  the  spectrum  shape  estimating  procedure  used  herein  (i.e.,  the 
numerical  curve-fitting  spline  function  and  spectrum  iteration  routine),  it  can  be 
expected  that  essentially  the  same  result  would  be  obtained  with  the  use  of  any 
other  spectrum  interpolation  schemes  which  can  accurately  account  for  curvature 
in  the  spectrum  shape. 

Once  this  general  spectrum  shape,  corresponding  to  the  initial  levels 
measured  with  real  filters,  is  estimated,  corrections  or  adjustments  for  different 
filters  (i.e.,  perfect  filters)  or  different  weather  or  propagation  distances  can  then 
be  applied  to  the  spectrum  at  all  frequencies  and  the  corrected  spectrum 
integrated  again,  for  each  filter  band,  to  obtain  new  adjusted  or  corrected  band 
levels.  Based  on  the  techniques  outlined  in  this  report,  these  corrected  band  levels, 
and,  for  aircraft  signals,  corresponding  values  of  perceived  or  effective  perceived 
noise  level  can  be  expressed  in  terms  of  values  that  would  have  been  measured  with 
an  ideal  filter  or  with  a  standardized  version  of  a  practical  filter. 

The  report  also  outlines  how  one  can  estimate,  within  certain  constraints, 
corrections  to  a  measured  spectrum  to  account  for  the  influence  of  additive 
background  noise  in  combination  with  the  effects  of  finite  filter  slopes  and 
extreme  signal  spectrum  slopes. 

Following  this  introduction,  Section  2  defines  the  type  and  magnitude  of 
anomalies  that  occur  in  the  spectral  analysis  of  aircraft  noise  due  to  the  individual 
and  combined  effects  of  large  negative  slopes  in  the  measured  spectrum  and  finite 
slopes  of  the  filter  skirts  with  or  without  the  presence  of  background  noise  and 
discusses  alternative  procedures  which  have  or  might  be  used  for  correcting 
spectra  for  filter  and  spectrum  slope  effects. 

Section  3  defines  the  general  iterative  approach  used  to  estimate  the  true 
spectrum  shape  and  outlines  the  design  of  a  computer  algorithm  based  on  this 
approach  which  can  be  used  to  correct  measured  band  levels  for  the  effects  of 
finite  filter  skirts  and  extreme  spectrum  slopes.  As  a  demonstration  of  this 
technique,  it  is  applied  to  a  real  aircraft  spectrum  measured  at  one  distance  and 
weather  condition  and  the  one-third  octave  band  levels  from  this  same  source  that 
would  be  measured  with  the  same  analyzer  are  predicted  for  a  new  distance  and 
weather  condition.  This  new  band  level  spectrum  is  then  used  as  an  input  to  the 
iterative  algorithm  and  the  original  measured  band  levels  are  obtained  very  closely. 
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Section  4  presents  the  results  of  applying  this  approach  to  a  variety  of 
idealized  and  actual  aircraft  spectra  to  further  demonstrate  the  relative  accuracy 
of  the  technique  and  to  illustrate  the  estimated  differences  between  PNL  or  EPNL 
values  obtained  when  this  technique  is  used  instead  of  current  methods  for 
specifying  band  levels  and  band  level  corrections.  The  former  allows  one  to 
simulate  and  evaluate,  explicitly,  errors  due  to  finite  slopes  of  analysis  filters  and 
signal  spectra;  such  errors  are  only  implicit  in  current  methods. 

It  is  shown  that  under  typical  conditions  specified  in  the  report,  these 
differences  in  PNL  or  EPNL  values  are  usually  less  than  I  dB.  However,  when  high 
frequency  one-third  octave  band  spectrum  measurements  are  made,  or  corrected 
to,  large  distances  under  conditions  defined  roughly  by  the  simple  rule  of  thumb, 
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band  levels  can  be  in  error  by  much  more  than  I  dB. 

Finally,  Section  5  summarizes  suggested  procedures  and  algorithms  for 
correcting  measured  aircraft  noise  spectra  for  the  influence  of  finite  filter  and 
extreme  signal  spectrum  slopes. 

Appendices  follow,  which  contain: 

A.  Simplified  model  for  output  of  a  finite  filter  for  both  signal  and 
background  input  noise  components. 

B.  Supporting  details  pertinent  to  this  report  in  the  form  of  a  reprint, 
with  errata,  from  a  recently  published  paper  in  the  Journal  of  the 
Acoustic  Society  of  America. 

C.  Analysis  of  power  transmission  through  a  practical  fractional-octave 
band  filter. 


2.  ANOMALIES  IN  BAND  LEVELS  FROM  SPECTRAL  ANALYSIS  OF 
BROADBAND  NOISE 

Three  possible  ways  to  define  band  levels,  identified  below  and  illustrated 
in  Figure  I,  might  be  used  to  interpret  the  spectral  content  of  a  broadband  noise 
which  is  evaluated  with  a  typical  spectrum  analyzer. 

o  The  reference  band  level,  L^,  for  a  reference  white  noise  signal 
which  hos  a  constant  spectrum  level  equal  to  that  of  the  signal  at  the 
band  center  frequency  f 

o  The  ideal  band  level,  Lg|,  of  the  true  spectrum  as  measured  with  a 
perfect  filter  over  the  nominal  filter  bandwidth  f  |  to  f2* 

o  The  measured  band  level,  Lq,  of  the  true  spectrum  as  measured  with 
a  practical  filter  over  the  significant  power  transmission  of  the  filter 
which  is  in  the  frequency  range  f^  to  fj. 

The  first  form,  L^,  is  identified  here  since  it  is  often  implied  as  the 
effective  band  level  when  carrying  out  frequency  variable  corrections  to  signal 
spectra.  For  example,  one  way  to  express  and  compute  the  attenuation  between 
two  bands  of  noise  due  to  any  attenuation  process  whose  magnitude  vories  with 
frequency,  implies  this  format  for  expressing  band  levels  when  the  attenuation  can 
be  accurately  defined  only  at  single  frequencies,  such  as  at  the  band  center 
frequencies. 

The  second  form,  LBI’  would  presumably  be  the  desired  form  for  universal 
application.  It  is,  in  fact,  closely  approximated  in  many  situations  when  the  signal 
spectrum  slope  is  not  large  and  spectrum  analysis  filters  are  employed  which  have 
very  steep  skirts. 

The  last  form,  Lg,  represents  the  actual  band  level  that  is  measured  with 
real  filters. 

Each  of  these  three  forms  can  be  described,  mathematically,  in  the 
following  fashion. 

LBo  =  Ls(fc)+  |0  log  [f2-f|]  ,  dB  (I) 


4 


Spectrum  Level,  dB 


Frequency  (log  scale) 


a)  White  Noise 


Figure  I.  Conceptual  Illustration  of  Three  Alternative  Ways  to  Interpret  a 
Band  Level  Centered  at  Frequency  f  with  a  Spectrum  Level  L  at 
this  Frequency  (*cross-hatched  area  equals  indicated  level) 
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(2) 


where 


spectrum  level  at  frequency  f 
Ls(f)  *  I0log10[p2(f)/p2] 

PQ  =  normal  reference  pressure  =  20  jxPa 

f|,  fc,  f2  =  lower  band  edge,  geometric  center,  and  ipper  band  edge 
frequency  of  an  ideal  filter 
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TL(f)  =  power  transmission  loss  of  real  filter 

f|_>  fy  =  lower  and  upper  frequency  limits  of  effective  transmission 
range  of  filter  *  f  j/5  and  5  •  fy  respectively 


This  report  will  consider  differences  between  these  three  forms  for 
expressing  band  levels,  first  for  single  and  then  multiple  bands  constituting  a 
complete  spectrum.  Finally,  however,  it  will  consider  differences  in  the  single¬ 
number,  frequency-weighted  and  time-integrated  metric  for  aircraft  noise, 
Perceived  Noise  Level  (PNL)  and  Effective  Perceived  Noise  Level  (EPNL)  respec¬ 
tively  when  these  are  evaluated  with  either  of  the  last  two  forms  above  for 
expressing  band  levels,  namely,  the  ideal  and  practical  band  levels  LBI  ar>d  *~B’ 

In  developing  this  information,  it  will  become  clear  that  the  first  form  for 
expressing  band  levels  is  also  very  useful,  conceptually,  as  an  intermediate 
transition  medium  for  estimating  changes  in  band  levels  due  to  any  frequency- 
sensitive  attenuation  process.  One  such  process  is  obviously  atmospheric  absorp¬ 
tion,  which  plays  such  a  strong  role  in  shaping  the  time-varying  spectrum  of 
aircraft  noise  signatures  during  a  flyover. 
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Define,  now,  two  necessary  and  sufficient  "error"  terms  which  will 
constitute  the  building  blocks  upon  which  the  principal  results  in  this  report  will  be 
developed.  (For  now,  the  compounding  effects  or  errors  in  spectrum  analysis  cbe 
to  background  noise  are  ignored.)  Let  the  reader  not  dispair  at  this  point.  These 
"error"  terms  are  really  used  only  to  aid  in  understanding  the  basis  for  the  results 
obtained.  These  results  will  be  expressed  in  more  practical  ways  which  do  not 
require  keeping  track  of  these  error  terms  explicitly. 

o  Spectrum  Slope  Error,  A 

This  will  be  defined  as  simply  equal  to  the  difference  between  the  ideal  and 
reference  band  levels,  Lgj,  and  L^,  so  that 

AS  =  LBI  "  LBo  =  LBI  “  [LS(fc)o+  10  ,ogIO  (f2  "  V]  ,dB  (4a) 

or,  solving  for  the  spectrum  level,  Lg(fQ) 

Lg(f Q)  =  Lgj  -  10  log |Q  (f 2  -  f  |)  -  Ag  ,  dB  (4b) 

As  will  be  clear  later,  it  is  this  second  form  in  eq.  (4b)  which  is  more  useful 
in  developing  a  direct  and  unified  approach  to  account  for  finite  filter  and 
signal  spectrum  slopes.  Strictly  speaking,  this  quantity,  A  5,  is  not  really 
an  error  in  spectrum  analysis  but  is  so  designated  here  for  simplicity  and 
consistency  in  terminology. 

o  Filter  Error,  A  p 

This  can  be  considered  a  real  error  in  band  levels  and  is  defined  as  the 
difference  between  band  levels  measured  with  a  real  (Lg)  and  ideal  (Lg|) 
filter,  or 

A  F  =  LB  '  LBI  ’  dB 

That  is,  Ap,  is  the  difference  between  the  band  level  measured  with 
practical  spectrum  analysis  filters  having  finite  slopes  of  their  filter  skirts 
and  the  ideal  band  level  that  would  be  measured  with  a  perfect  filter  with 
zero  transmission  outside  of  its  nominal  passband.  This  error  is  much  less 
for  one-third  octave  than  for  full  octave  band  filters  and  is  most 
significant  at  high  frequencies  and  large  distances  from  a  source  when  the 
spectrum  slope  can  be  very  steep. 


The  next  section  will  outline  one  unified  approach  to  aircraft  spectral 
analysis  which  has  the  potential  ability  to  occount  for  both  of  these  errors 
simultaneously.  First,  however,  it  is  desirable  to  obtain  a  quantitative  picture  of 
the  conditions  under  which  these  two  types  of  errors  occur  and  their  magnitude 
under  these  conditions. 

It  will  be  shown  that,  in  general,  differences  in  weighted  overall  aircraft 
noise  levels,  expressed  in  terms  of  PNL  or  EPNL,  will  usually  be  small  -  seldom 
exceeding  I  to  2  dB  and  usually  much  less.  However,  when  aircraft  noise 
signatures  are  being  analyzed  for  detailed  evaluation  of  source  contributions,  large 
errors  can  occur  in  specific  one-third  octave  band  levels  that  can  be  quite 
important  and  deserve  more  careful  consideration. 

One  example  of  how  these  "error"  terms  might  be  applied,  conceptual/,  is 
shown  in  Figure  2.  This  illustrates,  schematically,  one  basic  aspect  of  aircraft 
flyover  spectrum  evaluation  -  the  prediction  of  how  a  given  measured  spectrum  at 
one  point  -  identified  in  Figure  2  as  the  "source"  -  changes  as  the  sound  propagates 
to  another  -  the  receiver.  This  common  source-path-receiver  calculation  of 
propagation  loss  (and  corresponding  spectrum  shaping)  is  applied,  twice,  when 
correcting  data  measured  under  one  condition  of  weather  and  propagation  path 
length  to  another  condition  with  different  weather  and/or  path  length.  The  only 
difference  is  that  the  calculation  starts  at  the  receiver,  works  backward  to  the 
presumed  source  under  test  conditions  and  then  back  to  the  reference  receiver  for 
reference  conditions. 

As  indicated  near  the  bottom  of  Figure  2,  it  is  current  engineering  practice 
to  compute  a  receiver  band  level  Lg(R)  directly  from  a  source  band  level  Lg(s)  by 
subtracting  a  single  band  attenuation  *LB  due  to  atmospheric  attenuation  which 
depends  only  on  the  weather  conditions  and  total  path  length.^  (Spreading  loss  is 
ignored  here,  since  it  is  not  frequency-sensitive.)  The  alternative  approach  taken 
in  this  report  makes  it  possible  to  account  indirectly  for  the  two  analysis  errors  ip 
arK*  AS’  ^'rs*  a*  ^  source  aad  then  at  the  receiver,  as  illustrated  in  the  figure. 
The  problem,  illustrated  in  Figure  2,  is  treated  in  detail  in  reference  14,  included 
herein  as  Appendix  B,  but  with  one  variation.  The  difference  between  4  5  at  the 
source  and  receiver  was  used  to  define  a  single  correction  factor  relating  the 
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Figure  2.  Illustration  of  the  Conceptual  Application  of  Filter  and  Spectrum 
Slope  Errors  Ap  and  Aj  Respectively  to  Estimate  the  Propagation 
Loss  Between  a  Measured  Source  Band  Level  Lg(S)  and  an  Estimated 
Receiver  Band  Level  Lg(R).  In  each  case,  the  measured  band  level, 
Lq  is  corrected  by  Ap  to  obtain  the  ideal  band  level,  Lgj,  assuming  a 
perfect  filter.  Then  the  spectrum  slope  error  A  5  is  applied  to  Lgj, 
along  with  the  bandwidth  term  10  log  Ap,  to  obtain  the  spectrum 
level  at  the  band  center  frequency.  The  difference  between  the 
spectrum  levels  at  the  source  and  receiver  is,  ignoring  spreading  loss, 
simply  the  atmospheric  attenuation  ®(fc)  •  /at  the  center  frequency 
f  over  the  path  length  /where  ®(fc)  is  the  atmospheric  absorption 
coefficient  at  this  frequency. 
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difference  between  absorption  loss  for  ideal  bands  of  noise  and  the  single  frequency 
value  0f( f  )  -  A  at  the  band  center  frequency. 

It  must  be  emphasized  that  the  preceding  discussion  of  methods  for 
computing  absorption  loss  for  bands  of  noise  is  only  provided  here  as  an  illustration 
of  the  conceptual  application  of  the  two  "error"  terms  Ap  and  Ag.  Nevertheless, 
it  is  appropriate  to  bring  out  the  aspect  of  atmospheric  absorption  loss  which  has 
such  a  strong  spectrum  shaping  influence  on  aircraft  noise  signatures.  Specifically, 
the  type  of  high  frequency  roll-off  in  aircraft  noise  propagated  over  long  distances 
is,  in  fact,  the  principal  basis  for  considering  potential  errors  in  aircraft  single¬ 
event  noise  metrics,  such  as  PNL  or  EPNL,  due  to  the  combined  effect  of  these 
filter  and  spectrum  slope  effects.  Consider,  now,  more  specific  quantitative 
aspects  of  these  error  terms. 

2. 1  Filter  Error  Due  to  Finite  Filter  Skirts 

The  well-known  filter  error,  Ap,  is  inherent  in  the  fundamental  process  of 
signal  spectrum  analysis  and  has  been  evaluated  extensively  by  others  including,  for 
example,  Beranek  and  Sepmeyer.  ’  A  detailed  mathematical  derivation  of 
this  inherent  error  in  filter  band  levels  is  presented  in  Section  III  of  Appendix  B. 
The  approach  is  based  on  the  detailed  treatment  in  Appendix  C  of  the  power 
transmission  loss  TL(f)  term  for  the  practical  filter,  included  in  eq.  (3).  This  power 
transmission  is  simulated  by  defining,  mathematically,  a  suitable  expression  for  the 
filter  response  TL(f)  which  closely  approximates  the  response  of  real  filters.  Then, 
by  numerical  integration  of  the  power  transmitted  through  this  filter  and  an  ideal 
rectangular  filter  with  the  same  band  center  frequency,  the  filter  correction  factor 
A  j  is  computed  from  equations  (2),  (3),  and  (5). 

Examples  of  this  inherent  error  in  spectral  analysis  will  be  illustrated  later 
with  hypothetical  and  actual  aircraft  spectra  to  show  that  it  can  produce  sub¬ 
stantial  errors  in  specific  one-third  octave  band  levels,  even  close  to  a  source, 
when  the  spectral  shape  is  quite  uneven.  Clearly,  any  data  evaluation  procedure 
which  makes  it  possible  to  account  for  this  error  directly  could  be  very  useful. 
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As  pointed  out  earlier,  this  more  obscure  "error"  is,  in  effect,  implicit 
when  one  applies  corrections  valid  only  at  single  frequencies  to  band  level  data. 
Recall  that  this  error,  defined  by  equations  (I),  (2),  and  (4a),  is  simply  the 
difference  between  the  band  level  of  a  spectrum  measured  with  a  perfect  filter 
(Lgj)  and  the  band  level  (L^)  that  one  would  obtain  if  the  spectral  density  over  the 
entire  band  were  the  same  as  at  the  band  center  frequency. 

Based  on  the  mathematical  details  developed  in  Appendix  B  for  the 
integrated  power  transmitted  through  an  ideal  filter  band  (see  equations  (4)  to  (6b), 
(10)  and  (I  I)  on  page  B-3  of  Appendix  B),  one  can  derive  the  following  expression 
for  A  j  for  the  case  of  an  input  signal  with  a  constant  band  level  slope  applied  to  an 
ideal  filter. 

AS=  l0l°9|0  [r^(r(m+l)/2-r’(m+,)/2)/(m+l)(r‘l)]  ’dB  (6a) 


where 


r  =  ratio  of  upper  to  lower  band  edge  frequencies  f2  and 
fj,  respectively,  for  an  ideal  rectangular  filter 

m+l  =  S/10  log  2 


S  =  band  level  slope  of  the  true  input  spectrum  in  dB/octave 
For  the  case  of  m  =  -I  or  S  =  0  dB/octave  (pink  noise), 
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2.3  Filter  and  Slope  Errors  for  Constant  Slope  Source  Spectra 

Utilizing  the  concepts  just  defined  in  the  preceding  two  sections,  it  is  desirable  to 
examine  these  two  spectrum  analysis  error  terms  for  the  case  of  a  signal  input 
representing  a  source  signal  with  a  constant  band  level  slope,  before  spectrum 
shaping  by  atmospheric  attenuation  has  occurred. 


n 


The  results  are  shown  in  Figure  3.  Part  (a)  shows  the  slope  error  Ag(S),  as 
defined  by  eg.  (4a),  at  the  source  for  full,  one-third  and  I /18th  octave  band  filters. 
The  latter  represent  the  effective  spectral  filtering  bandwidth  to  be  employed 
later.  (Note  the  scale  change  for  this  case.) 

Part  (b)  shows  the  filter  error,  as  defined  by  eq.  (5)  for  the  full  and  one- 
third  octave  band  filters.  Note  that  this  error  is  of  the  order  of  1/2  of  A^  for  the 
range  of  slopes  considered. 

For  most  aircraft  noise  signals  measured  close  to  the  source,  say,  within  75 
m  (~ 250  ft),  the  band  level  slope  is  usually  less  than  9  dB/octave  unless  pure  tone 
components  are  present.  In  their  absence,  it  is  clear,  from  Figure  3,  that  both  Ap 
and  Aj  would  be  very  small,  for  one-third  octave  band  filters.  Under  such 
conditions,  the  actual  "as  measured"  band  levels  could  be  considered  free  of  side¬ 
band  effects  due  to  the  filter  skirts  and  the  spectral  density  at  the  band  center  can 
be  closely  approximated  by  the  "white  noise"  approximation  implicit  in  eq.  (I).  As 
we  will  see,  this  ideal  situation  does  not  hold  when  significant  tones  are  added  to 
the  broadband  signal. 

2.4  Filter  and  Slope  Errors  at  a  Distant  Receiver 

Consider,  now,  the  change  in  Ap  and  A^  at  a  receiver  due  to  spectral 
shaping  by  atmospheric  absorption.  A  variety  of  constant  slope  source  spectra, 
with  the  addition  of  a  pure  tone  component  in  one  case,  have  been  selected  for  this 
evaluation.  These  idealized  source  spectra  are  illustrated  in  Figure  4.  The  curves 
show  source  spectra  with  slopes  ranging  from  +3  dB/octave  to  -36  dB/octave  -  a 
range  sufficient  to  cover  normal  source  spectra  plus  special  cases  involving,  for 
example,  pseudotones,  where  local  spectral  slopes  at  a  source  signal  may  be  quite 
high. 

It  was  assumed  that  the  source  spectra  were  measured  with  an  ideal  filter 
so  that  the  initial  spectral  content  at  the  signal  could  be  determined  exactly. 
Then,  by  propagating  these  known  spectra,  through  known  atmospheric  absorption 
losses,  one  can  redefine  a  new  attenuated  spectra  at  a  receiver  and  then  apply  the 
band  integration  concept  described  at  the  beginning  of  this  section  and  developed 
in  detail  in  Appendix  C.  In  this  case,  it  was  possible  to  compute  new  band  levels  at 
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Figure  4.  Idealized  Source  Spectra  with  Various  Band  Level  Slopes,  Including  One  with 
o  Superimposed  10  dB  Pure  Tone  Component,  Used  to  Evaluate  Filter  and 
Spectrum  Slope  Errors 
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the  receiver  assuming  either  an  ideal  or  real  filter  (with  finite  sidebands)  and  hence 
compute  the  filter  error  Ap(R)  at  the  receiver.  Similarly,  the  received  spectrum 
level  can  be  defined  exactly  so  the  slope  error  Ag(R),  at  the  receiver,  can  also  be 
computed.  In  both  cases,  these  new  values  are  computed  from  the  same  equations 
(4a)  and  (5). 

The  results  of  these  computations  are  shown  in  Figures  5  and  6.  Figure  5 
shows  Ap(R)  as  a  function  of  frequency  for  the  source  spectra  of  Figure  4  and  for 
two  weather  conditions.  The  first  is  a  standard  day  of  25°C  and  70  percent 
relative  humidity,  for  which  atmospheric  absorption  losses  tend  to  be  near  a 
minimum  within  the  FAA  Part  36  test  window.^’  ^  The  second  is  at  I5°C  and  35 
percent  relative  humidity  -  a  condition  for  which  atmospheric  absorption  losses 
tend  to  be  a  maximum  within  the  limits  on  temperature  and  humidity  allowed  for 
aircraft  noise  certification.^’ 

The  general  trend  in  4  p(R),  shown  in  Figure  5,  can  be  explained  by  three 
interacting  influences. 

o  A  p  increases  at  all  frequencies,  uniformly,  as  the  source  spectrum 
slope  increases,  as  expected  from  Figure  3. 

o  A  p  increases  at  high  frequencies  due  to  the  further  increase  in  signal 
slope  due  to  atmospheric  absorption  over  the  fixed  600  m  propagation 
distance  assumed  -  this  effect  is  clearly  dominant  in  Figure  5b. 

o  A  p  is  modified,  in  a  complex  way,  for  any  bands  close  to  or  including 
a  band  with  a  strong  tone  component. 

In  Figure  6,  using  only  three  of  the  source  spectra  evaluated  in  Figure  5, 

the  variation  in  A<-(R)  at  a  receiver  has  also  been  evaluated.  Again,  the  same 

general  pattern  that  occurred  in  Figure  5  is  evident.  Note  that  these  error  terms 

are  not  shown  beyond  an  upper  bound  of  6  dB  intentionally.  They  generally  reach 

such  a  magnitude  only  when  the  total  signal  attenuation  due  to  atmospheric 

absorption  is  very  high  resulting  in  a  received  signal  near  or  below  practical  limits 

2 

of  measurement  in  the  presence  of  normal  background  noise. 

One  point  should  be  clear  by  now  -  the  real  error  term  A  p  and  the  pseudo 
error  Ag  behave  in  a  complex  manner  with  the  measured  signal  spectrum 
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Figure  5.  Voriotion  in  Filter  Error  5  at  Receiver  after  Propagation  Over  600m  with  Two 
Different  Weather  Conditions,  from  Source  with  Various  Initial  Slopes  Shown  in 
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Figure  6.  Variation  in  Slope  Error  a.  at  Receiver  After  Propagation  Over  600m,  with  Two 
Oifferent  Weather  Conditions,  from  Sources  with  the  Various  Slopes  Indicated 


characteristics.  This  complexity  is  we  1 1 -recognized,  of  course,  and  is  one  reason 
why  almost  no  attempt  is  made  to  account  for  these  error  terms  when  analyzing 
aircraft  noise  data.  One  exception  is  the  simple  procedure  inherent  in  the  aviation 
industry  standard  for  predicting  air  absorption*^  which  calls  for  using  essentially 
the  lower  band  edge  frequency  for  evaluating  all  absorption  losses  for  frequencies 
equal  to  or  greater  than  5000  Hz.  This  procedure  is  equivalent  to  trying  to 
recognize  the  inherent  slope  error  A<-  term  and  compute  atmospheric  absorption 
losses  at  a  frequency  below  the  band  center  which  should  be  more  "representative" 
of  changes  in  the  band  energy.  The  procedure  is  probably  in  the  right  direction  to 
aid  in  minimizing  errors  when  analyzing  attenuation  of  bands  of  noise.  However,  a 
more  general  approach  to  the  problems  of  accounting  for  filter  band  and  signal 
slope  errors  seems  desirable.  Such  an  approach  is  developed  in  the  next  section. 
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3. 


SPECTRUM  ITERATION  SCHEME 


The  preceding  analysis  of  the  filter  and  signal  slope  errors  Ap  and  WQs 
possible  because  the  initial  spectrum  levels  were  known.  In  general,  when 
measuring  an  aircraft  signal  at  some  distance  from  the  source,  one  really  doesn't 
know  the  true  received  spectrum  very  well  in  all  cases  when  current  techniques  of 
signal  analysis  are  used,  i.e.,  one-third  octave  band  analysis  with  filters  designed  to 
meet  applicable  industry  standards.  (See  discussion  at  the  beginning  of  Appendix  C 
for  additional  details  on  this  point.)  While  it  would  be  possible  to  use  narrow  band 
analyzers  to  evaluate  aircraft  flyover  signals  and  hence  determine  the  true 
spectrum  shape  within  each  one-third  octave  band,  such  a  process  would  probably 
be  impractical  for  universal  application  by  the  aviation  industry. 

Instead  of  requiring  actual  narrow  band  analysis  to  determine  the  true 
signal  spectrum,  it  seemed  desirable  to  develop  a  method  of  estimating  the 
apparent  true  spectrum  for  any  received  signal  such  that  by  integrating  this 
apparent  true  spectrum  over  the  full  frequency  response  of  the  filter,  one  would  be 
able  to  replicate  the  original  signal  "as  measured"  with  real  filters.  Once  this 
apparent  true  spectrum  is  available,  one  can  then  proceed  directly  to  carry  out 
desired  frequency-dependent  operations  on  the  original  signal  from  the  following 
perspective. 

o  The  apparent  true  spectrum  can  be  integrated  over  a  simulated  ideal 
filter  to  obtain  the  true  band  level  Lgj  free  of  filter  errors. 
Alternatively,  a  standard  "real  filter"  can  be  synthesized  to  compute 
band  levels  based  on  accepting  a  known  filter  error. 

o  There  is  no  longer  any  need  to  apply  corrections  for  the  filter  or 
signal  slope  errors  since  these  are  essentially  eliminated  once  the 
true  apparent  spectrum  is  known. 

o  The  apparent  true  spectrum  can  be  odjusted  for  off-reference  test 
conditions  of  measurement  distance  or  weather  using  atmospheric 
absorption  models  which  are  strictly  accurate  only  at  single  fre¬ 
quencies,  i.e.,  at  frequencies  corresponding  to  those  specified  for  the 
spectrum  levels. 
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A  spectrum  iteration  technique  to  accomplish  this  is  illustrated  conceptually  in 
Figure  7  where  it  is  shown  being  applied  to  estimate  the  true  source  spectrum  that 
would  be  measured  with  ideal  filters,  given  levels  measured  at  a  receiver  with  real 
filters. 

This  simplified  flow  chart  identifies  the  essential  steps  of  an  algorithm 
suitable  for  automation  in  a  computer  program,  details  of  which  are  described  in 
the  next  section.  Essentially,  the  problem  involves  estimating  the  true  spectrum 
levels  of  the  receiver  signal  and  iterating  this  estimate  until,  by  integration  of  this 
spectrum  (interpolated  at  any  desired  frequency  from  values  estimated  at  the  band 
center),  over  the  response  bandwidth  of  the  simulated  measurement  filter,  the  "as 
measured"  levels  can  be  duplicated,  within  a  definable  error  criteria.  Details  of 
this  process  are  covered  in  the  next  section. 

3.1 

Given  the  band  level  Lg(f.)  measured,  with  a  practical  filter,  at  each  i**1 
frequency  fj  in  the  range  f  j_  to  fy  we  wish  to  find  the  true  spectrum  level  function 
L  (f)  at  any  frequency  which,  when  integrated  over  the  effective  transmission 
response  range  of  each  mathematically  simulated  filter,  will  yield  the  same  band 
level  as  actually  measured.  This  actually  reduces  to  determining  a  spectrum  level 
L  (f :)  at  the  (geometric)  center  frequency  f.  of  each  band.  These  spectrum  levels, 
combined  with  a  suitable  interpolation  method,  are  then  used  to  estimate  spectrum 
levels  at  other  frequencies  and  hence  compute  band  levels  based  on  numerical 
integration  using  a  mathematical  representation  of  either  an  ideal  or  practical 
filter. 

It  should  be  recognized,  of  course,  that  such  a  procedure  for  estimating  the 
spectrum  levels  at  all  frequencies  within  a  given  band  is  subject  to  some  inherent 
limitations. 

o  The  estimated  spectrum  levels  can  be  uniquely  defined  for  each  band 
in  a  given  spectrum  with  the  use  of  any  one  suitable  iteration  and 
interpolation  scheme.  However,  any  change  in  either  the  iteration  or 
interpolation  scheme  may  produce  small  variations  in  the  estimated 
spectrum  levels. 
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Figure  7.  Spectrum  Iteration  Process  to  Estimate  True  Receiver  Spectrum  from  Measured  Receiver  Spectrum  and  thus  Allow 
Reconstruction  of  Source  Bond  Levels  Free  of  Errors  Due  to  Finite  Filter  Slopes  (Ap)  or  Spectrum  Slope  (A<j) 


o  The  estimated  spectrum  levels  are,  of  course,  a  smoothed  replica  of 
the  true  spectrum.  For  example,  they  represent  the  presence  of  a 
pure  tone  component  superimposed  on  a  broad  band  noise  by  means  of 
a  spectral  peak  with  a  finite  bandwidth  centered  near  the  actual  band 
center  frequency.  (An  example  of  such  a  representation  will  be 
shown  later.) 

Nevertheless,  the  iteration  and  interpolation  scheme  presented  in  this  report  has 
been  found  to  provide  meaningful  results  and  has  proven  itself  to  be  relatively 
insensitive  to  minor  variations  in  the  iteration  procedure,  thus  minimizing  the 
uncertainty  in  the  "uniqueness"  of  the  estimated  spectrum.  The  basic  steps  in  this 
procedure  can  be  defined  as  follows. 

3.1.1  Initial  Estimate  of  Spectrum  Level 

For  a  zero  order  approximation  to  the  spectrum  level  at  each  band-center 
frequency,  assume  that  the  noise  has  a  flat  spectrum  (i.e.,  constant  spectrum  level 
over  the  filter  band),  and  the  actual  measurement  filter  is  in  fact  ideal  (i.e.,  its 
transmission  is  unity  at  all  frequencies  within  its  nominal  passband  and  zero  at  all 
other  frequencies).  Thus,  this  zero-order  estimate  of  the  apparent  spectrum  level 
L  (f -),  or  simply  L  ,  at  the  band  center  frequency  will  be  (dropping  the  f.  index  for 

Lso“LB-IOI°9|o<4f)  P) 

the  nominal  filter  bandwidth 
[  2°  -  2”°  ]  f. 

’  1/2  for  full  octave  band  filters 
1/6  for  one-third  octave  band  filters 

nominal  band  center  frequency 
I()3n/I0  for  full  octave  band  filters 
Iqo/IO  fQr  octave  band  filters 

any  integer 
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3.1.2  Refined  Estimote  of  Spectrum  Levels 


Assume  three  spectrum  levels  or,  for  convenience  in  notation  at  this  point, 
their  equivalent  spectral  densities,  Sj_|,  S-,  and  Sj+|  at  frequencies  f._j,  fj,  and 
f.  respectively.  Further,  assume  the  approximate  value  of  the  band  level  is  the 
area  under  the  two  trapezoids  shown  in  Figure  8. 

Sj_—  Sj  between  f and  f ■ 

and  Sj+— -  Sj  between  fj+ and  fj 

where  f .+  and  f  .  are  the  upper  and  lower  bend  edges  of  a  constant  percentage  band 
filter  centered  at  f j. 

Further  assume  that  along  this  curve  of  spectral  density  vs  frequency,  the 
spectral  density  slope  changes  linearly  with  the  relative  frequency  so  that  at  any 
frequency  f,  the  spectral  density  is: 

,  .  A  '  W 

»f>  =  S|  (t  j  W) 

Letting  f  =  f._|,  then  between  fj_|  and  fj,  the  ratio  of  the  spectral  densities  is 


S(f 

-w 


S-fe) 


a-  +  b-(f -  i  /f  •) 
i  r i-r  r 


(9) 


or,  in  decibels, 


s<v  /V|\] 

6i-  =  10  '°9  5(Tj7JT  =  10  |°i  +  bi  (“f—j  I  '°9  T 
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Similarly,  between  fj  and  fj+|,  the  ratio  is 


S(fi+|)  /fl.1\ai  +  bl(Wfl 
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Spectrum  Level,  dB 


or,  in  decibels 
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Since,  for  constant  percentage  bands,  there  is  a  constant  ratio  r  given  by 


fi+l  fi 


=  r 


i-l 


equations  ( 1 0)  and  (II)  become 


=  10 

Oj  *  bj/r] 

log  r 

(12) 

=  10 

[Oj.bjf] 

log  r 

(13) 

Solving  for  Oj  and  b-(,  these  interpolation  parameters  for  the  i-th  band  are 


a.  = 


'  (r^  «'l)  10  log  r 


6.  -  t1  6. 

i+  i- 


(14) 


bi  = 


I 


(r  -  I)  10  log  r 


6.  .  -  6. 

i+l  i- 


(15) 


Applying  these  constants  to  the  general  expression  at  any  frequency,  the 
spectral  density  S(f),  or  the  equivalent  spectrum  level,  Lg(f)  =  10  log  [s(f)/P^jcan 
be  estimated  at  any  frequency  using  eq.  (8). 

This  process  is  used  to  estimate  the  spectrum  levels  at  the  filter  band 
edges  fj_  and  f.+  and,  as  indicated  in  Figure  8,  the  band  level,  assuming  a  perfect 
filter,  is  assumed  to  be  given  by  the  area  under  the  two  trapezoids.  (The 
integration  equations  for  the  power  in  each  trapezoid  have  already  been  cited  in 
Appendices  B  and  C.)  Essentially,  this  makes  it  possible  to  compute  the  slope  error 
A  5  for  this  spectrum  so  that  a  refined  estimate  of  spectrum  level  Lj(f.)  at  the 
center  frequency  of  the  i-th  band  can  be  obtained,  using  eq.  (4b).  But  this  refined 


25 


estimate  now  allows  a  second  estimate  of  the  interpolation  parameters  Oj  and  bj  to 
be  made  from  equations  (14)  and  (15)  using  new  values  for  the  spectrum  level 
differences  6.+  and  6._  between  adjacent  bands.  This  refinement  process  is 
repeated  twice  before  proceeding  with  the  next  step. 

3. 1 .3  Integration  of  Bond  Levels 

Now,  applying  the  same  band  integration  equations  cited  earlier,  the  total 
power  passed  by  a  simulation  of  the  real  filter  is  computed  numerically,  by  dividing 
each  band  into  b  constant  percentage  elements,  as  defined  in  detail  in  Appendix  C. 
A  value  of  b  =  6  will  be  shown  to  be  adequate. 

This  new  value  of  band  level  is  then  compared  to  the  original  "as  measured" 
value  and  the  error  compared  to  an  error  criteria.  (A  fixed  value  of  0.1  dB  was 
used  in  the  iteration  algorithm  applied  here.)  This  error  test  is  applied  to  each 
band,  one  at  a  time,  over  the  entire  frequency  range  and,  if  it  is  exceeded  for  any 
one  band,  the  entire  set  of  new  spectrum  levels  are  corrected  by  the  amount  of 
error  just  determined  and  the  entire  process  repeated,  as  indicated  in  Figure  7. 

Most  of  the  time,  two  or  three  iterations  were  required  for  real  aircraft 
spectra  to  achieve  stable  results.  The  outcome  of  this  process  is  then  a  set  of 
spectrum  levels  at  each  of  the  integration  elements  which  describe,  in  a  self- 
consistent  way,  the  actual  input  signal. 

One  minor  variation  was  allowed  in  this  process  for  some  cases.  If  the 
average  slope  over  a  given  band  was  less  than  a  predetermined  slope  criteria,  the 
band  integration  step  was  omitted  since  it  was  possible  to  compute  the  band  level 
directly  (i.e.,  Ap  was  expected  to  be  sufficiently  small  to  eliminate  the  need  for 
integration.) 

One  other  feature  explored  in  the  iteration  process  was  to  employ  an  error 
multiplier  when  correcting  the  old  spectrum  levels  by  the  error  in  band  levels. 
That  is,  the  new  spectrum  levels  at  the  band  center  frequency  were  defined  by 

LS(fc)new  =  LS(fc)old  '  C  (LB  new  ‘  LB  old) 
where  C  =  the  error  multiplier.  A  value  for  C  of  I  to  1.2  was  found  to  be 
satisfactory. 
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3.2 


Evaluation  of  Spectrum  Iteration 

A  test  spectrum,  shown  in  Figure  9a,  was  constructed  from  an 
actual  aircraft  flyover  signal  modified  to  contain  an  apparent  tone 
component  of  3150  Hz.  Figure  9b  shows  the  range  of  spectrum  levels, 
relative  to  final  base  case  values,  for  a  variety  of  minor  changes  to  the 
iteration  algorithm.  The  maximum  range  was  about  ±0.1  dB  -  consistent 
with  the  error  criteria  of  0.1  dB  for  the  iteration  process. 

Also  shown  in  Figure  9a  is  the  absolute  value  of  the  zero-order 
estimate  of  spectrum  level  for  each  band,  as  given  by  eq.  (7);  Figure  9b 
shows  the  deviation  of  this  estimate  from  the  base  case.  Note  that  the 
relative  error  in  this  value,  indicated  in  Figure  9b  by  the  diamonds  is 
much  greater  than  the  range  of  values  of  the  iterated  levels,  indicated 
by  the  shaded  areas. 

In  Table  I,  a  more  detailed  summary  is  shown  of  results  of  the 
minor  variations  in  the  iteration  algorithm.  The  adjustable  parameters 
were: 

o  The  number  of  bands  integrated  on  each  side  of  the  central 
i-th  band  (2  was  finally  selected,  thus  representing  each 
filter's  response  range  by  a  total  of  5  bandwidths). 

o  The  number  of  elemental  subdivisions  for  . integration  of 

each  band  (a  value  of  6  was  found  to  be  sufficient  for  good 
accuracy). 

o  The  error  multiplier  factor.  (A  value  of  1  was  finally 

selected. ) 

o  The  slope  criteria  used  to  short-cut  the  integration  of 

some  bands.  (This  was  usually  chosen  to  be  either  0  or  1 
dB,  and  caused  either  all  or  19  out  of  the  24  bands  in  the 
test  spectrum  in  Figure  9a  to  be  integrated  during 
iteration. ) 

The  results  of  these  various  cases  show  that: 

o  Three  iterations  were  usually  required, 

o  The  rms  error  in  the  iterated  band  levels  (re  the  "as 

measured"  values)  was  always  less  than  0.02  dB, 

o  Spectrum  levels  at  specific  frequencies  did  not  vary 

significantly  in  most  cases. 
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AL(f),  dB  L  (f),  Spectrum  Level,  dB  re  (20pPa) 


too  >000  O  KtOOO 

FREQUENCY  IN  HERTZ 


Figure  9.  Spectrum  Levels  Derived  from  Test  Case.  Part  a)  shows  iterated  values  (....)  and 

predicted  value  (o)  assuming  white  noise  approximation;  Port  b)  shows  range  (JZZ  777) 
of  L  (f),  relative  to  base  case  values,  for  variations  in  iteration  algorithm  and  the 
devotion  (o)  for  the  ideal  case,  shown  in  Port  a). 
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Table  1 

Summary  of  Spectrum  Iteration  Results  with  Test  Spectrum 
for  Various  iteration  Algorithm  Parameters 


Adjustable  Parameters 

Results 

Case 

No 

of1 

No 

.  of2 

Error3 

Slope^ 

'rL 

o 

• 

o 

Z 

RMS6 

Relative 
Spectrum  Level; 

dB 

Bands 

I  Elements  " 

Factor 

Criteria 

Iterations 

Error 

100 

1000 

4600 

Summed 

per 

Band 

Hz 

Hz 

Hz 

dB/Oct. 

dB 

A 

1 

t 

1 

3 

0.0092 

0.001 

0.002 

0.013 

B 

2 

A 

1 

3 

0.0148 

0.465 

-0.094 

-0.017 

C 

6 

0 

3 

0.0093 

-0.011 

-0.003 

0 

D 

1 

3 

0.0093 

0 

0 

0 

E 

.2 

1 

3 

0.0066 

-0.001 

-0.003 

-0.094 

F 

1 

2 

3 

0.0093 

N/A 

G 

.2 

2 

3 

0.0065 

-0.096 

-0.126 

-0.094 

H 

.5 

2 

2 

0.0337 

N/A 

I 

1 

8 

1 

1 

3 

0.0087 

N/A 

J 

» 

6 

1 

1 

3 

0.0094 

<0.001 

<0.001 

<0.001 

1 

2 

3 

4 

5 

6 
7 


No.  of  bonds  below  passband  included  in  integration  for  filter  response.  (Only  one  band  above 
passband  was  included  in  all  cases.) 

No.  of  integration  elements  into  which  each  band  is  divided. 

Error  in  Nth  band  estimate  multiplied  by  this  factor  to  correct  old  spectrum  level. 

Criteria  for  omitting  bands  from  integration  due  to  expected  low  filter  error  ftp. 

No.  of  times  iteration  algorithm  loop  is  exercised  to  achieve  required  accuracy. 

RMS  value  of  difference  between  calculated  and  actual  "as  measured"  band  levels  over  all  24  bands. 

Spectrum  level  at  three  frequencies  from  final  iterated  spectrum  taken  relative  to  Case  D  -  the  base 
case  used  for  remainder  of  study. 


As  a  final  test  of  the  concept,  the  apparent  spectrum  levels  developed 
were  used  to  compute  new  "measured"  band  levels,  simulating  a  real  filter  in  each 
case,  at  a  different  distance  and/or  weather  condition.  These  new  band  levels  were 
then  reinserted  back  into  the  spectrum  iteration  algorithm  and  used  to  compute 
their  own  spectrum  levels  and  thus  allow  correcting  these  new  "as  measured"  levels 
back  to  the  original  band  levels. 

Table  2  summarizes  the  results  of  this  exercise  showing  in: 

Col.  I  The  original  band  levels 

Col.  2  The  error  in  duplicating  these  by  iteration 

Cols.  3,  4  The  new  computed  band  levels  at  new  conditions 

Cols.  5,  7  The  error  in  duplicating  each  of  these  new  band  levels  when 

each  was  inserted,  in  turn,  back  into  the  iteration  algorithm 

Cols.  6,  8  The  final  resulting  "return  error"  or  difference  between  the 
original  and  recomputed  "as  measured"  levels 

Except  for  the  last  two  bands,  the  agreement  is  excellent.  This  residual 
error  is  attributed  to: 

1.  Lack  of  precise  knowledge  of  the  true  weather  conditions  upon  which 
the  original  test  data  were  based,  and 

2.  Difficulties  in  extrapolating  the  highest  frequencies  due  to  high 
spectrum  slopes. 

Nevertheless,  the  overall  results  are  considered  quite  satisfactory,  indi¬ 
cating  the  scheme  should  be  more  than  sufficient  to  test  for  errors,  due  to  filter 
sidebands,  in  PNL  or  EPNL  values  of  real  aircraft  spectra.  This  was  carried  out  in 
the  next  section. 
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Table  2 


Evaluation  of  Iteration  Algorithm  Starting  with  "Test”  Data  File, 
Generating  New  Spectra  at  Different  Distance  and  Weather  and  Using 
These  New  Files  to  Recompute  the  Original  Data  File 


II3SESBHMI 

1  ■■■■"■■  -  — .  —  ■  — ^ - - 

Test 

L _ 1250 

Column 

- 2 

3 

5 

6 

7 

8  — 

Temp.  °C 

o 

’■a- 

<M 

25 

o 

25 

24 

15 

24 

Humidity  % 

64 

70 

Bp 

70 

64 

35 

64 

Distance  m 

800 

600 

600 

750 

Frequency 

lb 

((f) 

lb 

lb 

*(f) 

e(f) 

c(f) 

(f) 

Hz 

L—i- _ _ _ _ _ _ 

- dD 

50 

79.9 

0.001 

79.91 

79.85 

0.01 

0.01 

0 

0 

63 

82.7 

0.001 

82.72 

82.63 

0.01 

0.01 

0 

0 

80 

85.5 

0.001 

85.53 

85.40 

0 

0 

0 

0 

100 

88.2 

0.001 

88.17 

87.98 

-0.08 

-0.15 

0 

-0.  16 

125 

89.1 

0.001 

89.10 

88.85 

-0.07 

-0.14 

0 

-0. 14 

160 

88.9 

0.001 

88.91 

88.59 

-0. 10 

-0.20 

0 

-0.20 

200 

88.3 

-0.002 

88.47 

88.09 

0 

0 

0 

-0.01 

250 

91.9 

0.002 

92.15 

91.76 

0.03 

0.02 

0 

0 

315 

91.8 

0.001 

92.  f  i 

91.76 

-0.05 

-0.10 

0 

-0. 10 

400 

91.8 

-0.002 

92.29 

92.08 

-0.02 

-0.03 

0 

-0.03 

500 

94.4 

0.001 

94.95 

94.95 

-0.10 

-0.21 

0 

-0.20 

630 

94.6 

0.001 

*5.33 

95.55 

-0.08 

-0.15 

0 

-0.  15 

800 

93.7 

-0.001 

94.60 

94.99 

-0.04 

-0.08 

0 

-0.08 

1,000 

94.1 

0.001 

95.06 

95.39 

-0.10 

-0.21 

0 

-0.21 

1,250 

92.7 

-0.001 

93.88 

93.74 

-0.05 

-0.08 

0 

-0.07 

1,600 

91.3 

0.001 

92.72 

91.49 

-0.09 

-0.09 

0 

0 

2,000 

88.5 

0.001 

90.27 

86.98 

-0.07 

-0.06 

0 

0 

2,500 

87.7 

0.002 

90.14 

82.73 

0.01 

0 

0 

-0.01 

3,150 

91.1 

0.012 

94.37 

81.86 

0.04 

0.05 

0 

0.01 

4,000 

81.5 

0.041 

85.82 

67.75 

0.07 

0.06 

0.01 

-0.02 

6,300 

72.2 

0.012 

83.06 

34.01 

-0.11 

-0.11 

0 

0.07 

8,000 

66.9 

-0.009 

83.72 

14.51 

-0.01 

-0.12 

0.32 

0.88 

10,000 

61.8 

0.001 

88.73 

-6.41 

0 

0.49 

0.61 

4.46 

RMS  Error,  dB 

0.006 

0.07 

0.15 

0.14 

0.93 

No.  of  Iterations 

3 

2 

6 

Col.  I  -  Hypothetical  "As  Measured"  Band  Levels  for  "Test"  Data  -  Real  Filter 


2  -  Difference  Between  Band  Levels,  Computed  from  Iterated  Spectrum, 

and  Col.  I. 

3,  4  -  Band  Levels  Computed  from  Iterated  Spectrum  at  New  Distance  and 

Weather  -  Real  Filter. 


5  -  Difference  Between  New  Band  Levels,  Computed  from  Iterated  Spec¬ 

trum,  Based  on  Treating  Col.  3  as  "As  Measured"  Levels  and  Col.  3. 

6  .  Difference  Between  Band  Levels  Computed  for  Original  Weather 

Conditions  of  Col.  1  and  Original  "As  Measured"  Levels  in  this  column 

7,  8  -  Repeat  of  Cols.  5  and  6  for  Different  Distance  and  Weather. 
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EVALUATION  OF  FILTER  EFFECTS  ON  PNL  OR  EPNL  VALUES 
Filter  Errors  in  PNL  Values 


Four  representative  aircraft  spectra  recorded  at  PNLTM  using  the  data 
collection  and  correction  procedures  described  in  Volume  I*  were  selected  for 
evaluation  of  any  PNL  errors.  (The  actual  data  were  corrected  very  slightly  to 
normalized  values  for  a  standard  day  and  a  measurement  distance  of  300  m.  This 
adjustment  was  generally  less  than  2  dB.)  These  are  shown  in  Figure  10.  Figure  1 1 
shows  the  high  frequency  portion  of  the  iterated  spectrum  level  for  one  case  (Data 
File  47  707  Approach)  compared  to  the  zero  order  estimate  of  spectrum  level  using 
eq.  (4b).  The  latter  appears  to  be  a  good  approximation  of  most  frequencies  except 
near  the  strong  tone  component  where,  for  example,  at  3150  Hz,  this  "white  noise" 
approximation  would  be  about  5  to  6  dB  higher  than  the  average  iterated  value  in 
this  band.  This  is  one  clear  example  of  how  the  zero  order  estimate  of  spectrum 
level  can  be  quite  misleading.  (The  filter  error  Ap  at  this  frequency  was  found  to 
be  3.4  dB.) 

In  Table  3,  a  detailed  listing  is  provided  of  the  "as  measured"  band  levels, 
Lg,  the  iteration  error  e  for  each  band  at  the  original  "as  measured"  distance  of 
300  m,  and  the  computed  filter  error  Ap  based  on  use  of  the  iterated  spectrum  in 
equations  2,  3  and  5. 

» 

In  addition,  the  values  of  new  computed  band  levels  Lgat  600  and  900  m 
are  given  along  with  the  corresponding  new  filter  error  Ap  and  the  difference 
A  between  this  new  band  level  and  a  value  that  would  have  been  computed  if 
only  the  single  band  attenuation  factors  specified  in  SAE  ARP  866A  had  been 
employed.  (For  the  sake  of  brevity,  only  the  higher  frequency  values  are  tabulated 
for  the  last  three  cases.) 

For  each  of  the  "as  measured"  and  computed  band  levels,  PNL  values  were 
computed  using  the  procedures  specified  in  FAR  Part  36,  Appendix  B,  augmented, 
where  necessary,  by  Reference  2 1 . 

These  PNL  values  are  summarized  in  Table  4  for  each  data  file  covering 
more  distances  than  shown  in  Table  3.  For  each  data  file,  the  computed  PNL  is 
listed  based  on  computations  with  the  iterated  spectrum  (whose  rms  error  is  also 
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Spectrum  Level,  dB  re  (20^ Pa)  /Hz 


75 


Figure  11.  Comparison  of  Initial  Estimate  of  Spectrum  Level,  Based  on  White  Noise  Approximation, 
Signified  by  Dashed  Line  and  Final  Estimate,  Signified  by  Solid  Line,  developed  from 
iteration  process  so  that  this  spectrum,  integrated  over  full  response  range  of  a  real 
filter,  produced  the  "as  measured  ".band  levels. 
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Table  3 

Iteration  Errors  and  Calculated  Filter  Errors  (Ap)  and  Error  in  Use  of  Simulated 
SAE  ARP  866A  Procedure  for  Four  Representative  Aircraft  Noise  Spectra 
Normalized  Initially  to  300  m,  25°C,  70  Percent  RH 

(o)  Data  File  35  -  727  T/O 


Distance,  300 

m 

Distance,  60 

Distance,  900  m 

Freq. 

«o<2> 

V3> 

V3) 

A  (5) 

asae 

A  (3) 

af 

a  (5) 

asae 

Hz 

dB 

dB 

dB 

50 

87.0 

-0.02 

0.98 

86.97 

-0.02 

-0.02 

86.95 

-0.03 

-0.03 

63 

98.2 

0.01 

-0.20 

98.19 

-0.20 

0.01 

98.17 

-0.20 

0.01 

80 

97.9 

-0.01 

0.16 

97.86 

0.14 

-0.01 

97.83 

0.16 

-0.01 

100 

101.9 

0 

-0.17 

101.85 

-0.17 

0 

101.81 

-0.17 

0 

125 

100.4 

0 

-0.02 

100.32 

r0.03 

-0.01 

100.25 

-0.03 

-0.01 

160 

99.9 

0 

0.03 

99.78 

0.03 

-0.01 

99.66 

0.02 

-0.01 

200 

103.1 

0 

-0.12 

102.92 

-0.12 

-0.01 

102.74 

-0.12 

-0.01 

250 

103.7 

0 

-0.10 

103.43 

-0.10 

0 

103.16 

-0.10 

0 

315 

102.7 

0 

-0.06 

102.29 

-0.07 

0 

101.89 

-0.07 

0 

400 

101.9 

0 

-0.08 

101.30 

-0.08 

0 

100.71 

-0.07 

0.01 

500 

101.2 

0 

-0.07 

100.34 

-0.08 

-0.01 

99.49 

-0.07 

0 

630 

100.5 

0 

-0.10 

99.33 

-0.10 

0.01 

87.17 

-0.09 

0.02 

800 

98.7 

0 

-0.03 

97.16 

-0.03 

0.01 

95.62 

-0.03 

0.02 

1,000 

98.4 

0 

-0.10 

96.45 

-0.09 

0.01 

94.50 

-0.09 

0.02 

1,250 

97.2 

0 

-0.08 

94.82 

-0.08 

0.01 

92.45 

-0.07 

0.04 

1,600 

95.4 

0 

-0.05 

92.55 

-0.05 

0.02 

89.71 

-0.05 

0.04 

2,000 

94.1 

0 

-0.07 

90.69 

-0.07 

0.01 

87.30 

-0.06 

0.04 

2,500 

92.8 

0 

-0.07 

88.67 

-0.07 

0.02 

84.57 

-0.05 

0.06 

3,150 

91.2 

0 

-0.08 

86.08 

-0.06 

0.06 

80.99 

-0.04 

0.15 

4,000 

88.4 

0 

-0.06 

81.83 

-0.01 

0.17 

75.34 

0.05 

0.42 

5,000 

84.5 

0 

0.01 

75.65 

0.10 

-1.04 

67.02 

0.25 

-1.76 

6,300 

82.0 

0 

-0.08 

69.56 

0.04 

-1.62 

57.61 

0.31 

-2.74 

8,000 

78.3 

-0.01 

0 

60.55 

0.30 

-2.20 

44.13 

1.26 

-3.07 

10,000 

76.3 

0 

-0.07 

49.94 

0.47 

-3.37 

27.23 

3.09 

-3.09 

PNL 

118.61 

-0.07 

115.91 

-0.07 

-0.07 

114.80 

-0.05 

-0.03 

Tcfcle  3  (Continued) 


(b)  Data  File  47  -  707  APP 


Distance,  300  m 

Distance,  120  m 

Distance,  600  m 

Freq. 

4(,) 

c(f)(2) 

A  <3) 

*F 

4<4) 

(3)  (5) 

*F  *SAE(> 

4W 

A  <3>  A  <5> 

F  °SAE 

Hz 

dB 

dB 

dB 

1,000 

79.9 

0 

-0.09 

81.09 

-0.01  0.01 

77.91 

-0.10  -0.03 

1,250 

81.4 

0 

-0.12 

82.84 

-0.11  0 

79.00 

-0.12  -0.01 

1,600 

81.9 

0 

0.29 

83.67 

0.32  -0.05 

78.96 

0.26  -0.07 

90.9 

0 

0.36 

93.06 

0.39  0.11 

87.31 

0.31  -0.17 

2,500 

96.8 

0.01 

-0.20 

99.27 

-0.19  -0.02 

92.72 

-0.20  0.07 

3,150 

83.2 

0.10 

3.38 

86.10 

3.21  -0.21 

78.68 

3.69  0.66 

4,000 

85.5 

0.01 

-0.03 

89.77 

0.01  0.23 

78.48 

-0.07  -0.28 

5,000 

88.2 

0 

-0.18 

93.71 

-0.17  0.82 

79.17 

-0.16  -1.22 

6,300 

83.1 

0 

0.17 

90.71 

0.13  1.12 

70.93 

0.43  -1.35 

8,000 

81.3 

0 

-0.12 

92.68 

-0.11  2.05 

63.30 

0.08  -2.40 

76.9 

0 

0.03 

93.69 

0  3.00 

51.47 

0.87  -2.40 

PNL 

114.37 

0.01 

0 

117.81 

0  0.32 

109.85 

0.03  -0.08 

(c)  Data  File  88  -  707  T/O 

Freq. 

Hz 

Distance,  300 

m 

Distance,  600  m 

Distance,  900  m 

1,000 

85.5 

0.08 

-0.07 

83.47 

-0.06  -0.07 

81.51 

-0.07  -0.07 

1,250 

85.5 

0.10 

-0.11 

83.02 

-0.10  -0.09 

80.63 

-0.10  -0.08 

1,600 

84.3 

0.06 

-0.07 

81.39 

-0.06  -0.04 

78.55 

-0.05  -0.02 

83.4 

0.15 

-0.08 

79.84 

-0.08  -0.15 

76.44 

-0.08  -0.12 

85.6 

-0.01 

1.32 

81.14 

1.14  -0.31 

76.71 

0.99  -0.60 

3,150 

96.0 

-0.02 

-0.20 

90.79 

-0.20  -0.03 

85.58 

-0.20  -0.06 

4,000 

87.5 

0 

0.91 

81.38 

1.14  0.62 

75.34 

1.40  1.32 

5,000 

78.9 

0.02 

0.39 

70.06 

0.60  -1.03 

61.68 

1.08  -1.60 

6,300 

81.1 

0 

-0.20 

68.45 

-0.16  -1.83 

56.15 

-0.08  -3.30 

8,000 

73.9 

0.04 

0.40 

57.05 

1.16  -1.30 

41.42 

2.47  -1.38 

10,000 

69.8 

-0.01 

-0.05 

44.00 

0.87  -2.81 

23.30 

5.44  -0.5 

PNL 

115.29 

0.01 

-0.03 

107.07 

-0.02  -0.08 

aiji 
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Table  3  (Concluded) 


(d)  Data  File  90  -  DC-10  T/O 


Distance,  300  m 

Distance, 

600  m 

Distance, 

900  m 

Freq. 

S0> 

c(f)(2) 

o 

■n 

w 

s<4) 

V3’ 

A  (5) 

asae 

■iw 

i> 

n 

to 

A  ® 

SAE 

Hz 

dB 

dB 

dB 

1,000 

83.7 

-0.09 

-0.08 

81.66 

-0.08 

-0.08 

79.71 

-0.08 

-0.07 

1,250 

82.7 

-0.07 

-0.06 

80.25 

-0.06 

-0.06 

77.87 

-0.05 

-0.04 

1,600 

82.7 

-0.08 

-0.07 

79.46 

-0.07 

-0.07 

76.60 

-0.07 

-0.07 

2,000 

82.6 

-0.10 

-0.09 

79.18 

-0.09 

-0.10 

75.66 

-0.10 

-0.10 

82.2 

-0.05 

-0.04 

77.98 

-0.05 

-0.07 

73.84 

-0.04 

-0.07 

3,150 

83.7 

-0.07 

-0.06 

78.37 

-0.08 

-0.15 

73.16 

-0.08 

-0.18 

4,000 

86.0 

0 

-0.07 

79.15 

-0.10 

-0.11 

72.40 

-0.10 

-0.12 

87.8 

-0.07 

-0.08 

78.58 

-0.11 

-1.41 

69.62 

-0.07 

-2.56 

6,300 

86.7 

0.01 

-0.11 

74.27 

-0.02 

-1.61 

62.25 

0.17 

-2.80 

8,000 

81.4 

-0,01 

0.07 

64,13 

0.55 

-1.73 

48.06 

1.59 

-2.24 

10,000 

75.5 

-0.02 

0.05 

50,45 

1.21 

-2.06 

29.73 

5.43 

0,21 

PNL 

110.85 

-0.05 

-0.08 

105.33 

- 0.05 

-0.25 

101.95 

-0.05 

-0.22 

Notes 

( 1 )  Lg 

1 

(2)  f(f) 

(3)  *F 

(4)  Lg' 
<5>*SAE 


Actual  "as  measured"  band  level  fed  into  spectrum  iteration 
routine. 

Difference  between  band  level  Lg1  computed  from  iterated 
spectrum  at  original  distance  of  300m  and  the  "as  measured" 
band  levels  at  this  distance. 

Filter  error  based  on  calculated  "as  measured"  band  levels 
with  real  filter. 

Band  levels  calculated  from  iterated  spectrum. 

Lg*minus  value  computed  from  calculated  "as  measured"  band 
level  using  only  simulation  of  SAE  ARP  866A  procedure. 
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Table  4 


Summary  of  PNL  Values  for  Four  Representative  Aircraft  PNLTM  Spectra 
PNL  for  Ideal  Filter  and  SAE  Method  Shown  Relative 
to  Calculated  “As  Measured"  Values  for  a  Real  Filter 


Data  File 

Temp. 

Humid. 

Filter 

Type 

rms(2) 

Error 

PNL  in  dB  at  Specified  Distance  in  Meters 

300(3) 

Meos'd 

75 

300 

450  600 

Calculated 

750 

900 

35 

727  T/O 

25° 

Real 

0.005 

118.61 

122.39 

118.61 

117.16 

115.91 

115.60 

114.80 

75* 

Idea) 

(3) 

•0.04 

0.07 

0.21 

0.07 

0.06 

o.os 

47 

707  APP 

Real 

0.020 

114.37 

118.82 

114.38 

112.00 

109.85 

N/A 

N/A 

Ideal.. 

(4) 

-1.03 

0 

•0.01 

-0.03 

. 

. 

(SAEJM 

•0.48 

-0.01 

0.06 

0.08 

- 

- 

88 

707  T/O 

Real 

0.011 

115.29 

119.33 

115.28 

112.98 

110.86 

108.91 

107.07 

Ideal 

(3) 

•0.051 

0.03 

0.03 

0.03 

0.02 

0.02 

(SAE) 

-0.31 

0.01 

0.07 

0.09 

0.09 

0.06 

90 

OC-IO  T/O 

Real 

0.017 

110.85 

116.98 

110.80 

107.77 

105.33 

103.20 

101.95 

Ideal 

(2) 

0.02 

0.08 

0.06 

0.05 

0.05 

O.OS 

(SAE) 

-1.14 

0.05 

0.25 

0.25 

0.2S 

0.22 

( 1 )  SAE  denotes  use  of  simulated  SAE  ARP  866 A  procedure  only  to  compute  basd  levels  at  other  distance  or  weather  conditions. 

(2)  RMS  value  of  difference  between  Iterated  end  actual  "as  measured"  band  levels  over  oil  24  bands  (no.  of  iterations 
in  parentheses). 

O)  PNL  far  the  actual  "as  measured"  spectrum  at  300  m. 
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listed)  assuming  a  real  filter.  Then,  the  value  that  would  have  been  measured  with 
an  ideal  filter  is  shown  relative  to  this  first  value.  Three  of  the  22  cases 
considered  show  an  effective  filter  error  in  PNL,  Ap  (PNL),  of  greater  than  0.2  dB. 
For  the  rest,  this  "error"  is  less  than  0.1  dB.  At  first,  this  is  surprising,  considering 
the  demonstrated  magnitude  of  the  filter  errors  Ap  in  individual  bands  for  these 
and  other  spectra  shown  earlier.  However,  the  explanation  seems  to  be  primarily 
that  an  individual  band  level  A  p  is  large  only  when  the  band  level  itself  is  already 
low  -  too  low  to  contribute  significantly  to  PNL. 

4.2  Filter  Errors  in  EPNL  Values 

Two  representative  fly-by  time  histories,  selected  from  the  data  defined  in 
Volume  ij  were  selected  for  processing  with  the  spectrum  iteration  algorithm, 
along  with  other  methods  for  fly-by  spectral  analysis  which  do  not  account  for 
filter  errors  but  do  attempt  to  account  for  missing  high  frequency  bands  due  to 
excessive  ambient  background  noise  by  various  techniques.  Results  of  this  trial  are 
summarized  in  Table  5.  Again,  it  is  clear  that  EPNL  values,  for  which  filter  errors 
were  eliminated,  listed  on  the  first  line  of  Table  5,  are  not  substantially  different 
from  the  values  computed  by  the  other  techniques.  Each  of  these  techniques  are 
defined  below. 

Reference  Method 

This  method  was  used  to  calculate  ideal  "Reference"  spectral  time 
histories.  The  procedure  was: 

1.  For  those  bands  within  5  dB  of  the  acoustic  ambient,  energy 
subtraction  of  the  ambient  from  the  signal  was  performed.  For  the 
approach  case  (flight  // 1 8),  a  total  of  5  bands  were  modified  using  this 
method.  For  //35,  none  were  modified. 

2.  For  those  bands  within  5  dB  of  an  electrical  noise  floor,  corrections 
were  made  using  linear  extrapolation  from  adjacent  frequency  bands 
(within  the  same  time  interval  that  the  ambient  violation  occurred). 
For  flight// 18,  two  bands  were  modified  using  this  method.  For 
flight  #35,  seven  bands  were  modified  using  this  method.  The 
frequency  extrqpolotion  was  performed  on  "normalized"  spectral  time 
histories.  These  consist  of  spectra  corrected,  for  inverse  square  law 
and  air  absorption  losses,*^  to  a  constant  propagation  distance  of 


Table  5 


Comparison  of  EPNL  Values  from  Two  Flights  Computed  by  Employing 
Spectrum  Iteration  Technique  and  Other  Conventional  Methods 


_ FPM 

' 

Method  Used 

Approach'  ' 
(#18) 

Takeoff(3) 4 

(#35) 

Reference 

With  Spectrum  Iteration^ 

104.19 

115.64 

Without  Spectrum  Iteration 

104.20 

114.90 

FAR  Part  36 

Within  5  dB(l) 2 

104.18 

114.99 

Within  2  dB 

104.21 

105.01 

Baseline 

(Wyle#l)i2) 

Up  to  4  Bands/Interval 

104.25 

115.03 

Up  to  8  Bands/Interval 

104.29 

115.10 

Refined 

Up  to  4  Bands/Interval 

104.24 

114.89 

(Wyle  #2) 

Up  to  8  Bands/Interval 

104.14 

114.87 

(1)  Margin  between  brckground  noise  and  signal  for  application  of  band  elimina¬ 
tion  per  Volume  II. 

(2)  EPNL  values  computed  as  defined  in  Volume  I. ' 

(3)  Measured  fly-by  time  histories  identified  in  Volume  ij 

(4)  The  EPNL  values  computed  with  this  method  were  based  on  "as  measured" 
spectra  corrected,  for  finite  filter  bandwidth  and  sideband  errors,  by  the 
spectrum  iteration  method  described  in  Section  3  and  corrected  for  ambient 
level  violations  by  the  method  defined  in  the  text. 
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60  m  (197  ft)  for  each  l/2-second  time  sample.  The  "normalized" 
spectra,  including  the  extrapolated  band  levels,  were  then  "unnorm¬ 
alized"  back  to  the  "as  measured"  distance  before  computing  EPNL 
values. 

Actuolly,  two  versions  of  this  Reference  method  were  used.  The  first 
method,  identified  on  the  first  line,  removed  the  filter  errors  by  employing  the 
spectrum  iteration  technique  defined  in  Section  3  on  the  "as  measured"  spectra 
before  making  the  additional  corrections  for  ambient  level.  The  second  Reference 
method  simply  omitted  the  spectrum  iteration  correction  so  filter  errors  remain. 

FAR  Part  36 

This  method  specified  that  all  bands  within  5  dB  or  2  dB  (two  different 
runs)  were  to  be  deleted  from  the  PNL  calculation  and  hence  EPNL  calculation. 

Baseline  Method  (Wyle  Method  //  I ) 

This  method  specified  that  all  (high  frequency)  bands  within  "X"  dB  of  the 
ambient  level  (electrical  ambient)  be  extrapolated  at  a  rate  of  -9  dB/octave.  "X" 
was  chosen  such  that  the  maximum  number  of  violating  bands  per  time  interval  did 
not  exceed  either  four  or,  for  the  second  run,  eight.  For  example,  for  the  approach 
case  (flight  //  18),  to  obtain  no  more  than  four  violations/time  interval,  a  signal-to¬ 
rtoise  ratio  of  16  dB  was  used  as  the  violation  criteria  Bands  affected  by  an 
acoustic  ambient  were  not  modified. 

Refined  Method  (Wyle  Method  //2) 

This  method  was  a  more  complex  form  of  the  "Baseline"  Method.  SPL 
values  were  first  "normalized"  to  a  distance  of  60  meters  for  each  I /2-second 
sample  of  the  spectral  time  history,  for  standard  day  weather  conditions,  and  then 
missing  bands  filled  in  by  extrapolation,  at  a  rate  of  -6  dB/octave  for  the  takeoff 
case  and  0  cB/octave  for  the  approach  case.  Acoustic  ambient  violators  were 
unaffected.  The  "corrected,  normalized"  levels  were  then  "unnormalized"  by 
converting  band  levels  back  to  "as  measured"  distances  in  order  to  compute  PNL 
and  EPNL  values  for  standard  day  conditions. 

4.3  Evaluation  of  Background  Noise  Effects 

Further  details  on  the  potential  compounding  effect  of  background  noise  on 
filter  errors  is  discussed  in  detail  in  Appendix  A  for  idealized  spectra  and  the 
reader  is  referred  to  that  section  for  details  and  conclusions. 
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5. 


CONCLUSIONS 


The  results  of  this  study  of  potential  errors  in  PNL  or  EPNt  values 
due  to  finite  energy  passed  by  filter  skirts  has  demonstrated  the  following: 

1.  For  the  test  cases  evaluated,  most  showed  a  very  small  error 
(less  than  0.1  dB)  due  to  finite  filter  effects. 

2.  However,  a  few  cases  showed  an  error  in  PNL  of  about  1.1  dB. 

3.  The  spectrum  iteration  scheme  developed  herein  has  proven  to  be 
a  valuable  tool  for  evaluating  filter  effects  and  is  a  promising 
means  of  eliminating  or  minimizing  the  substantial  filter  errors 
in  individual  bands  which  were  documented  in  many  of  the  cases 
considered. 

Thus,  as  an  improved  tool  for  critical  diagnosis  of  aircraft  noise  signatures, 
this  scheme  may  have  utility  in  the  industry.  Further  refinements  in  details 
of  the  algorithm  may  be  desirable  before  it  should  be  universally  adopted. 
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APPENDIX  A 


SIMPLIFIED  MODEL  FOR  SPECTRAL  ANALYSIS  OF  BROADBAND  NOISE 
CONSIDERING  BOTH  SIGNAL  AND  BACKGROUND  NOISE  INPUT  COMPONENTS 

I.  Introduction 

The  following  simplified  analysis  is  carried  out  to  provide  an  overall  view  of 
trends  in  the  measurement  of  aircraft  noise  when  effects  of  both  finite  filter 
slopes  and  background  noise  are  considered.  As  defined  later,  and  as  indicated  in 
Figure  A- lt  the  spectral  densities  of  the  signal  and  noise  elements  are  all  assumed 
to  have  slopes  defined  by  constant  exponent  relationships. 


Frequency,  Hz  (Log  Scale) 


Figure  A-l.  Signal  and  Background  Noise  Spectral  Density 


The  signal  and  two  background  noise  elements  are  introduced  into  the 
simplified  model  of  a  spectral  analyzer  system  in  the  way  illustrated  in  Figure  A-2. 
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Figure  A-2.  Conceptual  Diagram  of  Sound  Measurement  System  Illustrating 
Introduction  of  Signal,  S(f),  Pre-filter  Background  Noise,  N|(f), 
and  Post-Filter  Background  Noise,  ^(f),  into  Measurement  Chain 

All  elements  of  the  system  ere  assumed  to  have  unity  gain  over  the  nominal 
pass  band  of  the  spectrum  analysis  filter.  The  noise  and  signal  components  are  also 
assumed  to  be  uncorrelated.  Thus,  the  signal  power  observed  at  the  indicator 
(assuming  the  post-detection  noise  floor  of  the  system  is  well  below  the  measured 
signal  and  background  noise)  will  be 

P(f)=  [$(f)-N,(f)]  [T(f)  ♦N2(f)] 

It  should  be  noted  that  the  background  noise  component  defined,  in  Figure  A-2, 
as  post-filter  background  noise  N7(f)  is  essentially  the  same  as  specified  on  pages  4 
and  5  of  Volume  II  in  this  report  series,  as  "fixed  electric  background  noise"  (N^)  • 
Good  design  practice  in  any  spectrum  analyzer  would  dictate  the  need  to  strongly 
suppress  this  background  noise  component  that  is  added  to  the  signal  between  the 
filter  output  and  detector.  This  becomes  clear  when  one  recalls  that  the  total 
power  of  this  post-filter  noise  component  is  the  integral  of  its  noise  spectral 
density  over  the  entire  frequency  range  (i.e.,  greater  than  50  to  10,000  Hz)  of  the 
analyzer  while  the  power  in  the  signal  alone  is  roughly  equal  to  the  integral  of  its 
spectral  density  over  just  the  filter  bandwidth.  Hence,  to  avoid  contaminating  the 
filtered  signal  level,  the  spectral  density  of  the  post-filter  noise  (^(f)  in  Figure 
A- 1  or  A-2  must  be  substantially  lower  than  the  filtered  signal  spectral  density, 
$(f)  *T(f).  The  evaluation  in  this  c^jpendix  will  make  this  clear,  but  in  this  analysis 
of  background  noise  effects,  emphasis  is  placed  on  what  is  defined  here  as  pre¬ 
filter  background  noise. 
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Based  on  the  simple  model  defined  in  Figure  A-2,  the  following  defines  the 
frequency  variation  for  the  signal  and  background  noise  spectra  and  the  band  pass 
filter  response.  The  latter  is  a  very  simplified  version  of  the  frequency  response  of 
a  real  band  pass  filter;  this  simplification  was  chosen  for  convenience  in  the 
analysis  that  follows.  Therefore,  as  a  result  of  this  simplified  model  for  the  filter 
response,  we  will  be  primarily  concerned  with  relative  changes  in  the  measured 
band  level  due  to  the  introduction  of  pre-  or  post-filter  background  noise.  These 


relative  changes  in  output  are  considered  reasonable  estimates  of  the  true  change 

in  measured  level  when  both  background 

considered. 

noise 

and  finite 

filter 

slopes  are 

Let  S(f)  =  Spectral  Density  of  True  Signal 

= 

sc(f/fc)m 

(A- 1) 

Nj(f)  =  Spectral  Density  of  Pre-filter 
Background  Noise 

= 

N,(f/f/ 

N2(f)  =  Spectral  Density  of  Post-filter 
Background  Noise 

= 

N2(f/f/ 

1  * (f/f  L)q 

,f 

fL 

T(f)  =  Transmission  Response  of  Filter 

= 

1 

’fL 

f  fU 

Mf/fy)^ 

fU 

where  Sc,  Nj  and  N2  =  spectral  densities  of  the  signal  and  noise  elements  at 

the  geometric  center  frequency  fc  of  the  filter  band  in  question. 

A  more  convenient  form  for  the  spectral  density  slope  exponents  m,  n,  and  p 
and  the  filter  slope  exponent  +q  is  provided  by  the  corresponding  slopes  in  dB  per 
octave.  This  modified  form  can  be  defined  as  follows. 

If  the  signal  spectral  density  is  defined,  according  to  the  first  relationship 
above,  by  Sc(f/fc)m,  the  band  power,  B(f),  will  be 

B(f)  =  /S(f)df  =  Scfc(f/fc)m+l/(rfn+l)  (A-2) 

and  the  corresponding  band  level,  Lg,  will  be 
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(A-3) 


LB  =  10  lo9|0  [B(f)/Bref 
♦  IO(m+i)  logj0(f/fc) 


i0l°9|0  [scfc/(m+l)(Bref)] 


Thus,  for  any  octave  frequency  interval,  (f / f c)  =  2,  so  that  the  corresponding 
change,  m',  in  bond  level  over  this  frequency  interval  will  be 

m' =  iO(m+ 1 )  109^(2) '  ,  dB/octave 

and  the  corresponding  slope  exponent  quantity  (m+l)  becomes 

m+l  =  m'/IO  log  |  q(2)  (A-4) 

The  some  relationship  will  hold  tor  converting  the  band  level  slope,  in  dB/octave 
for  the  noise  elements  to  their  respective  slope  exponents,  n  and  p. 

A  similar  analysis  method  con  be  u^ed  to  define  the  slope  exponent  for  the 
filter  side-band  slopes  -  the  only  change  is  that  no  integration  is  involved  so  that, 
for  a  slope  of  q'  dB/octave  for  ft, ter  skirts,  the  equivalent  slope  exponent  q  will  be 

q  -  eq'./IO  log  (  q(2)  (A-5) 

For  typical  one-third  octave  band  filters,  q'  can  be  very  roughly  approximated  by  a 
slope  of  80  dB/octave. 

2.  Band  Power  Without  Background  Noise 

A.  Nominal  Passband  Power  for  Signal 

Consider  first  the  measured  band  levels  without  the  influence  of  background 
noise.  The  power  output  Bj  of  the  filter  between  the  nominal  filter  cutoff 
frequencies  fy  and  fy  for  a  signal  with  a  spectral  density  S(f)  proportional  to 
frequency  raised  to  a  constant  power  m  can  be  defined  by  the  integral: 


B,  =  J  S(f)df  Scf qJ  (f/fc)m  d(f/fc) 

fL  fL 

^c*c  '  *  (fi_/fc)m+*]//(m+l)  ,m£-l  (A-6) 

'  Scfc  ine(fU/fL) 
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Incorporating  the  relationships  between  fy,  fy  and  and,  for  convenience, 
expressing  the  total  power  relative  to  the  power  BQ  out  of  the  ideal  filter  with  a 


constant 

spectral  density,  S  input  at  the  band  center  frequency  fc,  eq.  A-l 

becomes 

Bl  =  Bo  Cl 

(A-7) 

where 

Bo  ■  Sc  [»U-*l]  =Sc  [<c  'h-,c  r'"] 

=  S  f  R 
c  c 

rr(fe)(m+l).r-Mm+l)l  /(m+l)R 

m  £-1 

C,  =  or 

,ne[r]/R 

m  =  -1 

and 

„  x  u  0  fc  -fc 

r  =  ^U'^L’  ^  =  r  ~ T- 

B. 

Lower  Sideband  Power  (By)  for  Signal 

Following  the  same  method,  the  power  (By)  in  the  lower  sideband  of  the 
filtered  signal  can  be  defined  as  follows.  (The  lower  sideband  is  assumed  to  be  just 
2  octaves  wide  extending  from  the  lower  band  edge  frequency  fy  to  fy/4.) 


By./  S(f)  T(f)  df  =  Scfc  /  (f/fc)m  (f/fc)q  (fc/fy)q  *f/fc)  (A-8) 
fy/4  fy/4 


where 

r('fe)  (m+ 1 }  [  I  -  <y*)m+q+ 1  ]  j  (m+q+  l)R  ,  m+q  k  - 1 

CL  = 

rq^  lne  ^4  j  /  R  ,  m+q  =  -I 

R  =  -  rh 
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C.  Upper  Sideband  Power  (By)  for  Signal 


Similarly,  for  the  upper  sideband,  the  signal  power  output  By  from  a  practical 
filter  is  assumed  to  be  represented  by 

/fU 

BU  =  Scfc  f  <f/fc)m  (f/g^  (fc/fu)"q  d(f/fc)  (A-9) 


rU 


BoCU 


where 


'U  = 


f(ft)  (m+l)  ^m-q+l  _  ,j  /[(m^|)Rj 
rq/2lne  [*]/R 


m-q  =  -l 


3.  Bond  Levels  of  Background  Noise 
A.  Pre-Filter  Noise 

For  the  pre-filter  background  noise  (i.e.,  background  noise  introduced  into  the 
system  before  the  spectrum  analysis  filter),  the  observed  band  level  can  be  defined 
in  precisely  the  same  way  as  for  the  signal.  Assuming  the  noise  spectral  density  is 
defined  by  N | (f )  =  Nj(f/fc)n,  then  equations  A-7  to  A-9  can  be  applied  directly  by 
simply  substituting  Nj  for  SQ  and  n  for  m.  Thus,  the  nominal  pass  band  power  B'j, 
lower  sideband  power  B1^,  and  upper  sideband  power  By  can  be  specified  by 

B'|  =  B'q  C'|  ,  power  of  pre-filter  noise  within  ideal  passband  (A-lOa) 

B'y  =  B'0  C'l  ,  power  within  lower  sideband  of  practical  filter  (A-lOb) 

B'y  =  B'q  O  ,  power  within  upper  sideband  of  practical  filter  (A-IOc) 

where 

B'o  *  NI«U  -  fL» 

and  C'|,  C'y,  C'y  are  the  same  as  defined  for  equations  A-2,  A-3,  and  A-4  with 
n  substituted  for  m.  Neglecting  the  special  cases  involving  n  =  I,  etc. 
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C*,  =  ^r(fc)(n*l).r-(fc)(n+l)j^  [(n+|)Rj  ,n*-l 

C'L  =r<-W<n+l>  [|.(fc)(n+<’+l)]/  [H+I)R]  ,  n+q  h  - 1 
C'y  sr(*)<n+,)  [4(^+,)-l]  /  [<n-q+l)R] 

B.  Post-Filter<  Bockqround  Noise 

The  background  noise  introduced  after  the  filter  is  present  in  the  measured 
output  over  the  entire  frequency  spectrum.  However,  it  will  be  assumed  that  for 
any  one  measured  band,  this  component  of  background  noise  is  present  over  2 
octaves  on  each  side  of  the  nominal  passband  of  the  one-third  octave  band  filter. 
Clearly,  this  implies  overlapping  in  the  accounting  for  this  noise  element  for 
adjacent  passbands.  However,  this  is  exactly  what  would  occur  since  each 
measured  band  level  is  considered  to  be  the  sum  of  all  signal  and  background  noise 
components  present.  Thus,  the  effective  band  power  of  this  noise  component,  when 
measuring  any  one  nominal  band,  will  also  be  broken  down  into  three  parts  as 
defined  below.  Again,  we  make  use  of  the  preceding  equations  A- 10  and  A- 1 1  by 
substituting  p  for  n  and  setting  q  =  o  to  eliminate  the  filter. 


B"l  ■ 

B"0  C", 

,  power  of  post-filter  noise  within  ideal  passband 

(A- 12a) 

b"l  = 

B”o  C”L 

,  power  within  frequency  range  of 
lower  sideband  of  practical  filter 

(A- 12b) 

B"u  = 

B'o  C"o 

,  power  within  frequency  range  of 
upper  sideband  of  practical  filter 

(A- 12c) 

where 

B"o 

=  N2(fu 

-fL) 

N2  =  Spectral  density  of  post-filter  noise  at  center  band  frequency  fc 
C",  =  [r(W(P+,)-r(-W<P+1)]/[(P+l)R]  ,p*-l 

C"L  ==  r(’W  (p* 1  *  [l  -  (y*)(p+ 1  )]/[(p+  DR]  ,  p  k  - 1 

C'y  =  r(fc)(p+l)  [4(p+l)- |]/[(p+DR]  ,  P  M 


(A- 1 1  a) 
(A- 1  lb) 
(A- lie) 
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4.  Combined  Siqnol  and  Background  Noise 

It  is  convenient  to  define  the  total  power  By,  measured  for  just  one  band,  in 
terms  of  the  following  matrix  of  band  powers  for  the  three  regimes  of  the  filter 
response  (i.e.,  lower  and  upper  sideband  and  ideal  passband  and  the  signal  and  noise 
elements  in  each  regime. 

ToTal  Measured  Power 

Lower  Ideal  Upper 


Sideband 

Passband 

Sideband 

Total 

Signal 

B  Ci 

0  L 

BoCl 

BaCU 

Bs 

Pre-filter  Noise 

B'  C\ 

0  L 

CD 

O  * 

O 

B’o  CU 

bn, 

Post-filter  Noise 

B"o  CL 

B"0  C", 

B’’o  C"U 

s 

Total 

bl 

Bi 

BU 

bt 

Only  the  outlined  cell,  Bq  Cj,  represents  the  power  in  the  true  signal  within 
the  ideal  or  nominal  passband  limits  of  the  filter.  All  the  other  elements  of  this 
matrix  represent  unwanted  added  power  segments  which,  under  certain  circum¬ 
stances,  can  cause  substantial  errors  to  occur  in  the  measurement  of  the  desired 
true  band  power.  Consider,  now,  some  quantitative  results  which  illustrate  such 
situations. 

A.  Results 

There  are  five  parameters  involved  in  portraying  the  results  of  this  appendix  - 
the  slope  exponents  m',  n’  and  p'  for  the  signal  and  noise  elements  respectively  and 
the  relative  levels  N|  and  N2  of  the  rxsise  elements.  Only  a  few  of  the  many 
possible  cases  are  considered  mere  in  detail.  A  summary  table  provides  a  brief 
overview  of  the  trends  over  a  wide  range  of  all  of  these  five  variables. 

Figure  A-3  illustrates,  for  six  specific,  cases,  the  relative  levels  in  each  of  the 
cells  in  the  preceding  matrix  of  signal  and  noise  elements  which  make  up  the  total 
measured  level.  The  cases  illustrated  cover  combinations  of  signal  and  background 
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Figure  A-3.  Illustration  of  Computed  Level  for  Each  Element  of  Total 
Measured  Output  from  a  One-Third  Octave  Band  Analyzer  Based  on 
the  Simplified  Model  Defined  in  the  Text. 


The  output  components  due  to  the  signal  I  ,  J  ,the  pre-filter 
noise  Y  /  /  A  ,  ond  the  post-filter  noise  POOP  are  given  for 
several  values  of  their  respective  slope  and  level  parameters.  The 
total  measured  level  is  signified  by  the  dashed  line  on  top  of  each 
figure.  The  total  difference  e  between  this  level  and  that  due  to 
the  signal  only  within  the  ideal  filter  passband  is  the  sum  of  the 
designated  components  ec?  due  to  the  added  signal  output  from  the 
filter  sidebands  and  eK,  due  to  the  background  noise. 


the  background  noise. 


noise  spectrum  slope  parameters  of  m'  =  0  and  -48  dB/octave,  n'  =  0  and  *6 
dB/ octave,  and  p'  =  0  and  +6  dB/octave,  respectively.  Values  of  the  noise  level 
parameters  are  N|  =0,-10  and  -20  dB  and  N2  =  -20  and  -40  dB.  The  center  bar 
represents  the  contribution  to  the  measured  level  within  the  ideal  filter  passband 
for  the  signal  and  both  noise  elements  (i.e.,  the  pre-  and  post-filter  noise),  while 
the  bars  on  either  side  present  the  corresponding  outputs  of  the  lower  and  upper 
filter  sidebands  (or  the  additional  +2  octave  range  included  in  the  post-filter 
background  noise).  The  reference  or  zero  level  for  each  of  these  elements  of  the 
total  measured  output  is  the  output  level  of  an  ideal  filter  for  a  white  noise  input. 
The  total  measured  level  is  designated  by  the  dashed  line  at  the  top. 

For  example,  for  case  a,  with  m1  =  -48  dB/octave,  the  portion  of  the  measured 
output  contributed  by  the  signal  within  the  lower  filter  sideband  is  about  4  dB 
above  the  reference  level;  the  corresponding  passband  level  for  the  signal  is  2.2  dli 
above  the  reference  level,  and  the  signal  output  from  the  upper  sideband  is  -18  dB 
below  the  reference  level.  (Again,  it  is  emphasized  that  these  values  are  expected 
to  have  the  correct  relative  magnitude  but  are  not  necessarily  representative  of 
actual  absolute  values  due  to  the  simple  model  chosen  for  the  filter  response.)  In 
this  case,  the  large  negative  slope  of  the  signal  input  spectrum  causes  the  lower 
sideband  output  of  the  filter  to  exceed  the  actual  passband  signal  output  by  nearly 
2  dB.  The  corresponding  filter  outputs  for  the  pre-  and  post-filter  noise,  indicated 
by  the  single-  and  double-shaded  bars,  are  symmetrical  about  the  passband  outputs, 
due  to  their  zero  band  level  slope,  and  the  levels  are  well  below  the  corresponding 
signal  output  as  expected  for  their  relative  spectrum  levels  of  Nj  =  -20  and  N2 
=  -40  dB  respectively. 

The  difference  e  between  the  total  measured  level  (signified  by  the  dashed 
line)  and  the  ideal  signal  band  level  within  the  filter  passband  is  the  sum  of  the 
designated  error  components,  e^,  due  to  the  presence  of  the  background  noise  and 
ep,  due  to  the  added  signal  output  from  the  filter  sidebands.  While  this  latter  error 
component  is  generally  larger  than  will  be  found  in  octuol  filters  (kje  to  the 
simplified  model  employed  for  the  filter  response,  it  is  only  shown  here  to  provide 
a  benchmark  for  evaluating  the  additional  error  due  to  the  presence  of  pre-and 
post-filter  background  noise. 


For  example,  consider  cases  (d)  and  (e),  which  differ  only  in  the  slope  of  the 
most-filter  background  noise  spectrum.  In  each  case,  the  spectrum  level  of  this 
noise  is  20  dB  below  that  of  the  signal  at  the  band  center  frequency  (i.e.,  ^  =  -20 
dB),  but  when  the  band  slope  of  this  noise  component  is  increased  from  0  to  +6 
dB/octave,  the  predicted  error  due  to  the  presence  of  background  nose  increases 
by  about  0.8  dB  to  1.54  dB.  Clearly,  the  requirement  that  post-filter  background 
noise  be  insignificant  suggests  that  its  spectrum  level  be  substantially  more  than  20 
dB  below  that  of  the  signal  level. 

Table  A-l  summarizes  the  values  of  the  error  e^,  as  defined  above,  for  values 
of  the  signal  slope  parameters  m'  and  n'  of  0,  -24  and  -48  dB/octave  and  0  and  -6 
dB/octave  respectively  as  well  as  for  values  of  the  pre-filter  noise  level  parameter 
N|  of  0,  -5,  -10,  and  -20  dB.  As  suggested  at  the  beginning  of  this  appendix,  good 
design  practice  for  a  spectrum  analyzer  would  dictate  that  the  post-filter  noise 
would  be  suppressed  sufficiently  to  be  entirely  negligible  so  it  has  been  ignored  in 
this  final  summary. 

The  principal  conclusions  that  can  be  drawn  for  this  appendix  are  as  follows: 

1.  Assuming  the  post-filter  background  noise  is  negligible,  the  error,  e^,  due 
to  the  remaining  background  noise,  is  defined  by  the  difference  in  energy 
between  the  total  measured  signal  and  the  background  noise  alone. 

2.  When  the  signal  and  pre-filter  spectrum  slopes  are  the  same  (i.e.,  both 
zero),  the  background  noise  error  e  ^  is  simply  equal  to  the  usual 
increment  in  level  when  two  uncorrelated  signals  are  added  on  an  energy 
basis. 

3.  The  background  noise  error  decreases  substantially  when  the  signal  has 
a  steep  negative  spectrum  slope  -  a  condition  that  will  occur  when 
analyzing  high  frequency  bands  of  aircraft  noise  measured  at  relatively 
long  distances  and/or  in  conditions  of  high  air  absorption. 

4.  To  avoid  significant  errors  in  the  measured  band  level  of  a  signal,  due  to 

background  noise  (i.e.,  <  0.1  dB),  the  spectrum  level  of  the  pre-filter 

background  noise  should  be  at  least  20  dB  below  that  of  the  signal  at  the 
band  center  frequency.  For  practical  purposes,  in  aircraft  noise  analysis, 


Calculated  Error,  e  ^  in  dB,  in  Total  Measured  Band  Level 

Due  to  Presence  of  Background  Noise  Introduced 
Prior  to  One-Third  Octave  Band  Anolyzer  Filter 


m',  Slope  of 
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Level 

dB/Octave 

Nj,  Noise 
Spectrum  Level* 
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n',  Slope  of  Noise 

Band  Level 
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-6 

0 

-48 

-20 

0.01 
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-24 

-20 
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0.04 
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-5 
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♦Spectrum  level  of  pre-filter  noise  at  center  frequency  of  band 
relative  to  signal  spectrum  level. 

♦♦Equal  to  normal  correction  for  energy  addition  of  two 
uncorrelated  signals. 


this  translates  to  a  predictable  requirement  that  the  aircraft  signal  to 
background  noise  ratio  should  be  >.  least  20  dB  in  each  one-third  octove 
band.  This  is  an  expected  result  which  certainly  did  not  require  the 
preceding  analysis  to  establish  it.  However,  when  this  signal-to-noise 
ratio  cannot  be  met,  as  is  often  the  case  in  analysis  of  high  frequency 
bands  for  aircraft  noise,  then  one  can  conclude  from  the  preceding  the 
following  significant  and  not  immediately  obvious  conclusions. 

5.  The  background  noise  error  is  nearly  independent  of  the  slope  of  the 

pre-filter  background  noise  spectrum  over  the  range  of  slopes  (0  to  +6 

dB/octave)  likely  to  be  encountered.  (As  shown  in  Volume  II  of  this  report 

2 

series,  the  slope  of  the  acoustical  background  noise  is  typically  approxi¬ 
mately  -4.5  dB/octave  while  the  spectrum  slope  of  the  pre-filter  back- 

2 

ground  noise  (or  electric  background  noise  in  the  parlance  of  Volume  II) 
can  vary  from  at  least  -3  to  +3  dB/octave.) 

6.  Thus,  the  total  band  level  of  the  signal,  including  the  power  in  the  filter 
sidebands  but  excluding  any  background  noise  contribution  can  be  deter¬ 
mined  by  simply  subtracting,  on  an  energy  basis,  the  pre-filter  (or  pre¬ 
detection)  background  noise  band  level  from  the  total  measured  band 
level.  This  noise-free  signal  band  level  can  now  be  used  with  the 
spectrum  iteration  procedure  outlined  in  Section  3  of  the  main  text  to 
determine  the  approximate  shape  of  the  true  spectral  density  of  the  signal 
thus  enabling  one  to  compute  true  signal  band  levels  and/or  adjust  the 
signal  levels  for  other  distance  or  weather  conditions  and  eliminate  the 
measurement  error  due  to  finite  filter  sidebands  in  the  presence  of 
moderate  pre-filter  background  noise. 

7.  To  insure  that  this  last  procedure  is  carried  out  correctly  for  all  types  of 
signals  contaminated  by  pre-filter  background  noise,  it  may  be  desirable 
to  actually  process  the  "as  measured"  signal  through  the  spectrum 
iteration  routine  first  in  order  to  accurately  define  the  true  "signal  ♦ 
noise  spectrum.  Then,  by  applying  an  energy  subtraction  of  the  true  noise 
spectrum,  one  obtains  a  true  noise-free  signal  spectrum  without  further 
iteration  which  is  ready  for  applying  frequency-sensitive  corrections 
and/or  computing  true  band  levels  simulating  ideal  filters.  (See  Section  3 
for  a  further  discussion  of  this  concept.) 


A-13 


APPENDIX  B 


INFLUENCE  Of  ATMOSPHERIC  ABSORPTION  ON  THE 
PROPAGATION  OF  BANDS  OF  NOISE 


By 

L.  C.  Sutherland 
H.  E.  Bass 

Reprint  of  Paper  Published  in 
Journal  of  Acoustical  Society  of  America 

September  1979 


ERRATA 


Page 

B-1  Abstract  -  Fourth  line  from  bottom;  change  "4  to  40"  to  read  "4  and  10" 


B-2 


B-6 


Section  1.  -  Third  line  after  Eq.(2);  change  to  read  "transmission  loss  of 
practical  filter  *“  f^/5  and  5  •  fj,  respectively  ..." 


Eq.  (26)  The  first  bracketed  term  should  read 


10  log 


nb 

Z  B 

LJ* 


rJ 


B-7  Figure  5  and  Figure  6.  In  titles,  substitute  Ap( S)  for  6^(S)  and  6p(R)  for  6j(R), 
respectively 

B-9  Section  V.  -  Fourth  line  from  bottom;  change  "distance"  to  "distant" 

B-10  Reference  16  -  Should  read  "American  Nationol  Standard  SI .  1 1— 1966(R— 1971 ).  . 


1 


Influence  of  atmospheric  absorption  on  the  propagation  of 
bands  of  noise 

L  C.  Sutherland 

Wfk  Laboratories.  El  Sepvndo,  California  9024S 

H.  E.  Bass 

Department  of  Physics  and  Astronomy.  The  University  of  Mississippi.  University.  Mississippi  38677 
(Received  1 1  January  1979;  accepted  for  publication  30  April  1979) 

At  a  band  of  tsoue  propagates  through  the  atmosphere,  atmospheric  absorption  attenuates  the  higher 
frequencies  in  the  band  more  rapidly  than  the  lower  frequencies.  For  large  propagation  distancea, 
atmospheric  absorption  shapes  the  spectrum  such  that  moat  of  the  acoustic  energy  is  in  Che  lower  part  of 
the  band.  When  this  occurs,  the  effective  atmospheric  absorption  coefficient  for  this  band  of  noise  is  less 
than  when  the  energy  is  evenly  distributed  over  the  bend.  When  i  real  filter  ie  used  to  examine  the 
received  spectrum,  the  problem  becomes  even  more  complex  since  the  filter  does  not  cut  off  completely  at 
the  nominal  lower-frequency  limit  for  the  band.  The  equations  governing  the  attenuation  and  spectrum 
analysts  of  noise  under  these  conditions  have  been  derived  and  a  numerical  integration  process  has  been 
used  to  obtain  quantitative  results.  It  is  shown  that  predicting  observed  values  of  atmospheric  absorption 
lam  liar  bands  of  noise  involves  two  correction  factors  which  are  added  to  single-frequency  abeorptioo  toes 
to  define  the  effective  attenuation  for  the  band  of  noise.  The  first  correction  factor  A,  accounts  for  the 
difference  between  the  attenuation  predicted  over  a  given  path  in  a  homogeneous  atmosphere  for  a  pure 
tone  at  the  geometric  center  frequency  of  an  ideal  filter  band,  and  the  actual  attenuation  over  the  tame 
path  for  the  same  ideal  band  of  noiee.  The  second  correction  (actor  L,  accounts  for  the  difference  between 
the  hand  level  measured  at  one  point  with  a  practical  filter  with  finite  transmission  outside  the  nominal 
filter  paes  band  and  the  true  band  level  measured  at  the  same  point  with  an  ideal  (rectangular)  filter. 
Representative  calculations  show  that  the  sum  of  these  corrections  become  substantial  when  the  product  of 
distance,  in  km,  and  the  square  of  frequency,  is  kHz,  exceeds  about  4  to  40  for  octave  and  one-third 
octave  band  filters,  respectively  Although  the  net  effect  of  such  corrections  will  generally  be  small  for 
overall  frequency  weighted-uotse  levels  (e  g.,  perceived  noise  level),  the  errors  in  di«p.r»in[  spectral  detail 
in  none  signal  urea  of  distant  sources,  such  as  aircraft,  can  be  important. 

FACS  numbers:  43.2*.Fp,  43.30  Vt 


INTRODUCTION 

Recent  experimental1*’  and  theoretical**’  work  have 
made  it  possible  to  compute,*  with  high  precision,  the 
acoustic  absorption  of  a  pure  tone  as  it  propagates 
through  a  still  atmosphere  free  from  the  influence  of 
a  ground  plane.*  Although  there  are  still  some  unre¬ 
solved  questions  concerning  the  contribution  to  the  at¬ 
tenuation  of  the  signal  by  turbulence,*  radiation,*  dust,10 
fog,"  or  finite  amplitude  effects,1*  the  major  contribu¬ 
tions  to  the  pure  tone  absorption  coefficient  are  now 
well  understood  and  documented.  In  practice,  however, 
one  is  seldom  confronted  with  the  propagation  of  pure 
tones.  Instead,  the  source  of  interest  typically  gener¬ 
ates  some  type  of  noise  spectrum  with  energy  spread 
continuously  over  a  range  of  frequencies.  It  is  the  ab¬ 
sorption  of  such  broadband  noise,  then,  rather  than  of 
a  pure  tone,  that  is  more  frequently  of  practical  inter¬ 
est,  particularly  for  the  analysis  of  air  to  ground  prop¬ 
agation  for  aircraft  noiae. 

Atmospheric  absorption  ia  an  increasing  function  of 
frequency.  As  a  result,  when  a  band  of  noise  propa¬ 
gates  through  the  atmosphere,  the  higher  frequencies  in 
the  band  are  attenuated  most  strongly.  If  the  spectrum 
is  Initially  flat,  after  propagating  some  distance,  it  will 
begin  to  exhibit  greater  and  greater  attenuation  at  high 
frequencies  as  either  distance  or  frequency  increases. 
This  spectrum -shaping  effect  of  atmospheric  absorp¬ 
tion  compounds  the  problem  of  spectral  measurements 


of  broadband  noise  signals  which  propagate  over  long 
distances  in  the  atmosphere. 

The  objective  of  this  paper  ia  to  outline  the  errors 
that  may  occur  and  outline  appropriate  correction  pro¬ 
cedures  when  one  evaluates  the  atmospheric  absorption 
loss  of  bands  of  noise  in  terms  of  the  pure  tone  absorp¬ 
tion  coefficient  <*(/,)  at  the  center  frequency  ft  of  a 
band.  As  outlined  in  Fig.  1,  two  basic  types  of  errors 
are  involved.  The  first,  designated  by  the  symbol  AJ(  is 
the  difference,  in  decibels,  over  a  fixed  path  at  length 
l  In  a  homogeneous  atmosphere,  between  the  absorption 
loss  lor  the  pure  tone  corresponding  to  the  center  fre¬ 
quency  cf  a  band  of  noise  and  the  loss  for  the  total  pow¬ 
er  within  the  same  band  of  noise  as  measured  with  an 
ideal  filter.  The  second  type  of  error,  designated  by 
the  general  symbol  A,,  equals  the  difference,  in  deci¬ 
bels,  between  the  total  power  passed  by  a  practical  fil¬ 
ter  and  the  power  passed  by  its  nominal  (ideal)  pass 
band.1*’1*  The  first  correction  factor,  A/,  will  vary  with 
the  nominal  filter  bandwidth  and  the  slope  at  the  source 
or  receiver  spectrum,  when  analyzing  propagation  lots 
from  a  source  to  a  receiver,  or  from  a  receiver  back 
to  a  source.  In  addition,  A,  will  vary  with  the  atmo¬ 
spheric  conditions  over  the  propagation  path  and,  unlike 
the  pure  tone  absorption  coefficient  a(fj  at  the  band 
center  frequency,  A,  will  also  vary  with  distance  along 
the  propagation  path.  The  second  correction  factor, 

A„  can  be  defined,  completely,  in  terms  of  only  the  fre¬ 
quency  response  at  the  filter  and  upon  the  spectrum 
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FIG.  1.  Illustration  of  the  application  of  two  con-eciiou  fac¬ 
tor*,  Ai  and  A^,  for  analyzing  propagation  loss  of  bands  of 
noise  A,  is  the  difference  between  the  pure  tone  loss 
(<»t/,MI  at  the  band  center  frequency/,  over  a  path  length 
I  between  a  source  IS)  and  receiver  Uf'.  and  the  to**  over 
the  same  path  for  an  Ideal  band  of  aoiae  (.entered  0,1  ih.  same 
frequency.  A^  is  the  difference  In  decibels  between  ti  e  power 
level  Lj  transmitted  by  an  actual  filter  with  finite  transmission 
outside  the  nominal  passband  and  the  power  level  lBl  -aosralt- 
ted  by  an  ideal  (rectangular  response)  filter  centered  no  'be 
tame  frequency. 


shape  of  its  input  signal.  However,  the  spectrum  shape 
of  the  filter  input  signal  at  a  receiver  is  also  intiuemied 
by  the  same  atmospheric  and  path  length  factor,  that 
influence  the  ideal  noise  band  to  pure  lone  correction 
factor  A,.  Thus,  in  many  cases,  it  is  convenient  to 
compute  one  composite  correction  term,  a,  which  var¬ 
ies  with  all  of  the  parameters  just  defined. 

The  fact  that  atmospheric  absorption  will  stupe  (he 
spectrum  within  a  band  of  noise  has  been  well  recog¬ 
nized  for  some  time.  For  example,  the  SAE  standard 
on  air  absorption”  requires  that  at  frequencies  above 
4  kHz,  the  pure  tone  absorption  coefficient  at  Die  lower- 
band  edge  frequency  should  be  used  to  compute  effec¬ 
tive  absorption  coefficients  for  octave  bands,  while  at 
lower  frequencies,  the  center  frequency  of  the  band 
should  be  used.  While  this  is  a  reasonable  rule  of 
thumb,  considering  typical  aircraft  noise  spectra,  it 
can  lead  to  significant  errors.1 

This  paper  will  present  a  mathematical  r  ,  reserva¬ 
tion  for  the  propagation  of  a  band  of  noire  through  the 
atmosphere  and  numerical  procedures  for  defining  the 
correction  factors  A,  and  A,  which  yield  quantitative 
results  not  dependent  on  approximations  such  se  used 
in  the  SAE  standard.  The  numerical  procedure  is  used 
to  compute  band  loss  coefficients  for  a  few  typical 
cases.  Due  to  the  many  variables  involved,  including 
weather,  initial  spectrum  shape,  type  of  filter,  and  dis¬ 
tance,  only  a  few  limited  cases  can  be  presented  »  ithin 
the  scope  of  this  paper.  Condensed  tables  are  also 
available  for  a  few  specific  cases.1'* 


t.  METHOD  OF  ANALYSIS 

In  geneial,  the  analytical  methods  used  to  define  the 
correction  factors  A,  and  A,  involve  the  computation, 
by  integration  over  a  specified  bandwidth,  of  the  rela¬ 
tive  level  (In  decibel s)  of  a  broadband  random  signal 
with  a  nonumlorni  spectrum.  The  analysis  does  not 
attempt  to  include  composite  signals  made  up  of  broad¬ 
band  noise  plus  one  or  more  pure  tone  components. 
However,  such  signals  can  be  handled  by  a  superposi¬ 
tion  process,  applying  the  methods  outlined  herein  for 
computing  the  absorption  loss  of  the  random  noise  por¬ 
tion  and  the  simpler  methods  applicable  for  computing 
absorption  losses  of  pure  tones.* 

The  specific  expression  which  defines  A,  can  be  stated 
as  follow  'throughout  this  analysis,  correction  factors 
are  derived  for  only  one  filter  band  but  they  can  obvi¬ 
ously  be  applied  to  any  number  of  bands  by  the  same 
techniques.) 

ftroyjd  *  lOlogJf  ip\(j)-\0-a'/",,xadf  j 

/  p\u) & ) ,  a) 

where  i  f.)  is  the  pure  tone  absorption  coefficient  at  the 
hand  center  frequency/,,  dB*m~‘;  l  is  the  propagation 
path  lengih.  m,  P\(f)  is  the  pressure  spectral  density, 
at  source  reference  position,  at  frequency/,  Pa'^Hz"1; 
and/,,  /,  are  the  nominal  lower  and  upper  cutoff  fre¬ 
quencies  oi  ideal  filter  band  centered  on  frequency/,. 

Hz. 


The  filter  correction  factor  A r{R)  at  the  receiver  is 

defined  t>v 

•SjriR)  ■-  10  log P\f)\0'rL'>'n(1  dl  j 

dB  ,  (2) 

where  Tlif)  is  the  power  transmission  loss  of  filter, 
dfi :  fr  .  t\:  are  the  effective  lower  and  upper  limits  of 
transmission  of  practical  filter  =»/,  8  and  8 •/,,  respec¬ 
tively,  Hz;  and  P‘r(f)  is  the  pressure  spectral  density  at 
input  to  receiver  filter,  Pa;,Hz"‘. 

A  comparable  expression  can  also  be  written  for  the 
filter  correction  factor  A r(S)  at  the  source  by  substitu¬ 
ting  the  subscript  s  for  rin  the  preceding  equation. 

Since  the  integrands  in  these  equations  cannot,  in  gen¬ 
eral.  be  expressed  in  a  simple,  integrable  closed  form 
valid  over  the  full  width  of,  say.  a  -J-odt3'-6  band  filter, 
i'  is  necessary  to  resort  to  a  process  of  numerical  in¬ 
tegration.  The  usual  case  of  interest  is  for  constant 
!>ei('en‘age  filters  for  which  the  ratio  of  the  upper  and 
iower  nominal  band  edge  frequencies/,  and/,  can  be 
specified  as  r -/,.'/,  and  the  geometric  center  frequency 
is  /,  -  (/,  jri)‘/1.  First,  the  bandwidth/,-/,  of  the  actual 
filter  is  divided  into  b  constant  percentage  rectangular 
filter  elements.  For  the ;th  element,  the  ratio  of  the 
upper  band  edge  frequency  f,.,  to  the  lower  band  edge 
frequency  /,  is  equal  to  rl/*  and  the  frequency /,  is  equal 
to 

/,-/,< Hz  »> 
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where  i  varies  from  1  at  the  lower  band  edge  at  the  en¬ 
tire  filter  to  (A ♦  1)  at  the  upper  band  edge.  Examine, 
bow,  the  application  of  this  numerical  integration  ap¬ 
proach. 

II.  DERIVATION  OF  A, 

Consider  first,  the  integral  for  the  power  in  an  ideal 
filter  band  at  the  source  as  specified  by  the  denomina¬ 
tor  of  Eq.  (1).  If  the  bandwidth /,.,  — /,  of  the  jth  ele¬ 
ment  within  this  ideal  filter  band  is  sufficiently  small, 
it  is  possible  to  assume  that,  over  this  narrow  frequen¬ 
cy  range,  the  pressure  spectral  density  at  the  source 
can  be  specified  by  a  simple  power  law  expression  of 
the  form 

////)"'.  Pa’-Hz-1,  «) 

where  P\jmP*,(/j),  the  pressure  spectral  density  at  the 
lower  edge  (/«/,)  at  the  jth  element,  and  equals  the 
exponent  which  defines  the  slope  of  the  spectral  density 
over  the  frequency  range /,-/,♦,.  This  simple  power 
law  is  equivalent  to  assuming  that  a  log-log  plot  cf  the 
spectral  density  versus  frequency  is  a  straight  line  with 
a  constant  slope,  over  the  frequency  interval  f,  to 
defined  by  the  exponent  mt.  The  total  power  B„  within 
this /th  elemental  band,  at  the  source  (corresponding  to 
the  mean  square  pressure  in  the  source  signal),  is  given 
by  the  integral  of  P\(f)  over  the  same  frequency  inter¬ 
val.  Thus,  with  Eq.  (4),  one  has 

)"'<?.  P»'-  15) 

With  Eq.  (3),  this  integral  for  the  total  power  in  the  >th 
filter  element  at  the  receiver  reduces  to 

Pa1,  n*,  *  - 1 ,  (ga) 
Pa1,  m, «  - 1 . 

Summing  over  the  entire  band  gives  the  total  band  pow¬ 
er  B,  at  the  source 


The  total  air  absorption  loss  over  the  source- recei¬ 
ver  path  l  can  also  be  given  to  a  good  approximation, 
over  the  same  narrow  frequency  range  /,  to /,.,,  by  a 
power  law  of  the  form, 

>**s,  (7) 

where  At  *A(/,)’'10""<0,‘,/,°,  the  fractional  loss  In 
power  by  absorption  at  the  lower  band  edge  frequency 
ft  at  the yth  filter  element,  and  n,  is  the  exponent  which 
defines,  approximately,  the  change  in  total  absorption 
with  frequency  over  the  narrow  frequency  range  of  this 
jth  element.  Then  the  pressure  spectral  density  at  the 
receiver  P\(f)  can  be  expressed,  with  Eqs.  (4)  and  (6) , 
as 

P*rV)  ~P\{f)'A(f)  •P\t(f//t)m>At{J //,)-!,  Pa*-  Hz-1 . 


Thus,  by  comparison  of  Eqs.  (4)  and  (S),  we  can  use 
Eq.  (6)  to  express  the  total  power  BrJ  (or  corresponding 
mean  square  pressure)  in  the  j th  elemental  band  at  the 


receiver  aa: 

'J>JfAj/<[ln.(r)]/b,  1. 

Again  summing  over  the  entire  band,  the  total  band  pow¬ 
er  at  the  receiver  is 


Br  •  B., . 


Each  of  the  elope  exponents  m,  and  n,  can  be  determined 
from  their  corresponding  log-log  relationships.  For 
example,  mj  Is  determined  by  setting  the  frequency  / 
«/iM  in  Eq.  (4),  taking  the  log  at  both  sides,  and  solving 
for  m,  to  give 

Substituting  the  spectrum  level  L,(f)  corresponding  to 
the  value  of  101ogloPi(/)  and  using  Eq.  (3),  this  reduces 
to 

mt  »  M,* i)  -  ^^)]/[UO/6)log10r] .  (11) 

Applying  the  same  process  for  nt,  using  Eqs.  (7)  and  (3), 

n,  •  [«</,♦.)  -<*(/,))•  l/[(\0/b)  log,„r] .  (12) 

The  values  at  a(f)  required  in  this  last  expression 
are  computed  with  the  methods  of  Ref.  6  at  the  two  ele¬ 
mental  filter  frequencies  ft  and /JM  and  for  the  average 
atmospheric  conditions  prevailing  over  the  path  length 
l.  (Note  that  this  path  can,  if  necessary,  be  just  one 
arbitrarily  small  incremental  segment  in  a  layered 
atmospheric  model.  A  different  correction  factor 
is  then  computed  for  each  such  elemental  path  segment 
and  for  each  frequency  band.) 

The  parameters  4  and  r  are  preassigned  constants. 
Based  on  the  numerical  integration  method  employed  to 
sum  powers  over  each  filter  band,  the  constant  b  (the 
number  of  integration  elements  per  band)  is  an  even  in¬ 
teger  and,  to  achieve  acceptable  accuracy  (integration 
error  <0.1  dB),  it  must  be  at  least  12  for  either  i  or  f 
octave  bands.  The  bandwidth  ratio  r  is  set  equal  to  109-1 
or  2  for  i  or-J  octave  band  filters,  respectively. 

All  that  remains  in  order  to  carry  out  the  integration 
over  all  b  segments  in  each  filter  band,  is  a  means  of 
specifying  the  spectral  density  or  corresponding  spec¬ 
trum  level  at  any  frequency  /  within  the  overall  filter 
bandwidth/,-/,.  Two  approaches  have  been  employed 
as  illustrated  in  Fig.  2.  The  first,  and  simplest,  was  to 
assume  that  the  slope  of  the  spectrum  level  curve  is 
constant  over  the  full  nominal  bandwidth  of  the  filter 
under  consideration.  The  second  approach  assumed  that 
the  spectrum  level  was  given  by  a  general  interpolation 
or  spline  function  C(f)  which  defined  a  spectrum  level 
curve  which  passes  exactly  through  the  true  spectrum 
levels  at  the  geometric  center  frequencies /,_,,  /,,  and 
/,.,,  of  the  »th  band  and  the  immediately  adjacent  bands. 

A.  Constant  slope  approach 

If  the  spectral  density  of  the  source  has  a  constant 
slope  over  the  entire  frequency  bandwidth  of  a  filter, 
then  it  can  be  defined  by  Eq.  (4)  letting  P*u  be  the  spec¬ 
tral  density  P\(fi)  at  the  lower  cutoff  frequency  (/,)  of 
the  filter  and  m,  becomes  simply  the  constant  slope  m. 
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FIG.  2.  Illustration  of  the  two  general  methods  considered 
to  defloe  the  spectrum  level  L,  If)  at  any  frequency  for  analy¬ 
sis  of  the  propagation  loss  for  the  >th  band  of  noise  with  center 
frequency  /,. 


A  more  convenient  definition  of  the  spectrum  shape  is 
provided  by  the  slope  (Sa)  of  the  filtered  band  levels  in 
dB  per  bandwidth.  Since  the  power  within  a  constant 
percentage  band  increases  directly  with  frequency  tu 
the  first  power  for  a  constant  input  spectral  density,  if 
the  spectral  density  varies  with  frequency  to  the  power 
m,  then  the  band  power  varies  with  frequency  to  the  pow¬ 
er  frn+1).  Thus,  if  the  ratio  of  upper  to  lower  cutoff 
frequencies  is  r,  then  the  slope  Sa  or  change  in  band 
level,  in  decibels,  over  a  frequency  interval  of  just  one 
bandwidth  (for  which  the  ratio  of  upper  to  lower  fre¬ 
quencies  is  r)  will  simply  be 

SB  =  101oglo(r)"*1  =  10 Ifn  *  1)  logl0r,  dB  •  (band)-1 . 

(13) 

Following  this  same  “constant  slope”  approach,  assume 
that  one  can  approximate  the  actual  frequency  variation 
in  the  atmospheric  absorption  over  one  entire  filter 
bandwidth  by  the  simple  power  law  given  by  Eq.  (7).  Re¬ 
placing  Aj  by  the  corresponding  absorption  loss  A ,  at  the 
lower  cutoff  frequency /,  of  the  filter,  then  at  any  fre¬ 
quency  within  the  band,  the  fractional  loss  due  to  ab¬ 
sorption  is 

•  (14) 

[This  is,  in  fact,  a  rather  crude  approximation  chosen 
here  for  simplicity  in  the  analysis.  A  better  approxima¬ 
tion  for  the  frequency  variation  in  atmospheric  absorp¬ 
tion  loss  over  the  path  l  is  equal  to  exp[-a,(/, /,)*•/], 
where  a,  Is  the  absorption  coefficient  (in  nepers  m'1 
instead  at  dB-m'1)  at  the  reference  frequency/,  and 
the  exponent  k  need  only  vary  from  0  to  2.]  Again,  a 
more  convenient  way  to  describe  the  variation  in  ab¬ 
sorption  as  a  function  of  frequency  by  this  approximate 
form  is  given  by  the  slope  SA  or  change,  in  decibels,  of 
the  absorption  loss  over  a  frequency  interval  of  one  fil¬ 
ter  bandwidth.  Thus,  the  parameter  SA  is  given  by 

SA  »  101og10(r)"**-10nlog10r,  dB- (band)'1 .  (15) 

From  Eq.  (1),  the  correction  factor  A,  can  be  expressed 

in  simpler  form  as 


A,  *a(/J'/tl01ogl0(Br/B,l,  dB,  (16) 

where  8r  and  B,  are  the  band  powers  at  the  receiver  and 
source  respectively  for  ideal  filters. 

Further,  by  transforming  the  first  term  into  decibel 
notation  and  employing  Eq.  (14),  it  can  be  expressed  as 

a(/c)-f--101oglo[A(/e)l*-101oglo[AiC/,//l)-J,  dB,  (17) 

where /,//,  =  rl/*,  the  ratio  of  the  geometric  mean  fre¬ 
quency  (/,)  to  the  lower  cutoff  frequency  (/,)  of  the  fil¬ 
ter. 

Now,  by  using  Eqs.  (6)  and  (9)  to  define  the  band  pow¬ 
ers  B,  and  B,  at  the  receiver  and  source,  respectively 
[dropping  the/  subscript  since  we  are  now  considering 
only  one  integration  element  over  the  entire  filter  band¬ 
width  (i.e. ,  b  -  1)],  and  including  Eq.  (17),  Eq.  (16)  be¬ 
comes,  for  the  nonsingular  case  (i.e.,m#-l  or  m-n 

A,  =  10  log  I0[B,/(B,A>r'-"1)] 

=  10 logl0[(r— 1 "♦*  -  l)(m  +  1)/ (r~n /a)(r",M  -  1  1)]  . 

(18) 

Employing  the  identities  r*-r~'  =  2sinh(ilnr)  and 
logl0r/lnr  =  log,0r,  Eq.  (18),  after  a  bit  of  algebra  and 
including  Eq.  (13)  and  (15),  becomes 

A]  =  10  loglo{Sa  sinh^S*  +Sa)]/(SX  +  Sa)*sinh(;3Sa)}  dB, 

(19a) 

where  0  =  1/20  log10r  =  0.11513. 

For  the  singular  cases  for  which  m  =  -1  or  m  -  n  =  -1, 
the  equivalent  for  Eq.  (19)  is  that  SB  =  0  or  SA*SB  =  0. 

For  these  cases,  it  can  be  shown  that,  for  Sa  =  0, 

A,  =  10  log10(sinh(/BJ  'RSj.  dB,  (19b) 

and,  for  SA+SB  =0, 

Ay  =  101og10[3Sa/sinh(3Sj)],  dB  .  (19c) 

Tlius,  simple,  though  approximate,  expressions  are  ob¬ 
tained  in  terms  of  just  two  variables.  SA  and  Sa,  for  the 
correction  factor  A,  which  specifies  the  difference  be¬ 
tween  the  single  frequency  absorption  loss  and  the  true 
value  for  an  ideal  band  centered  on  the  same  frequency. 
The  slope  SA  of  the  absorption  loss  can  be  estimated 
by  computing  the  values  for  the  total  absorption  loss 
at  frequencies  equal  to  one  filter  bandwidth  below  (/c/r) 
and  above  {Jc-r )  the  center  frequency  (/c)  of  the  filter 
in  question  and  computing  the  average  slope  by 

SA~[a(fc-r)-a(/'/r)]-l‘2<  dB-(band)*'.  (20) 

The  slope  of  the  source  spectrum  SB  can  be  approxi¬ 
mated  in  a  similar  fashion  or  estimated  by  graphical 
means. 

The  correction  factor  A,  for  an  ideal  [-octave  band 
filter,  as  defined  by  Eq.  (19),  is  presented  graphically 
in  Fig.  3  as  a  function  of  Sa  and  SA  in  the  form  of  con¬ 
tours  of  constant  values  of  Ar  While  the  values  of  A, 
shown  are  only  as  accurate  as  the  approximations  em¬ 
ployed  in  deriving  Eq.  (19),  they  provide  a  useful  first 
approximation.  There  are  several  characteristics  at 
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FIG.  3.  Cootour*  of  the 
value  of  approximate 
correction  factor  (A/).  - 


Fig.  3  of  interest.  First,  for  typical  conditions,  the 
correction  A,  is  not  large;  it  varies  from  about  6%  to 
8%  of  the  total  pure  tone  absorption  loss  a(/,)*l,  with 
the  highest  values  occurring  at  the  largest  distances. 
Second,  A,  is  usually  positive  since  for  most  full  scale 
noise  sources,  the  energy  in  a  band,  after  propagating 
through  the  atmosphere,  is  concentrated  at  frequencies 
near  the  lower  band  edge  where  atmospheric  absorp¬ 
tion  is  lower  than  at  the  band  center.  Thus,  the  true 
absorption  loss  for  the  band  is  less  than  the  pure  tone 
loss  of  the  band  center  frequency.  Note  that  for  this 
simplified  approach,  the  values  of  A,  in  Fig.  3  are  in¬ 
dependent  of  the  filter  bandwidth  since  the  slope  pa¬ 
rameters  SA  and  St  are  in  decibels  per  bandwidth.  In 
fact,  of  course,  the  actual  value  of  A,  lor  an  octave  band 
filter  will  be  substantially  greater  than  for  a  J-octave 
band  filter  for  the  same  actual  spectrum  and  absorption 
slopes  in  dB /octave.  We  can  conclude  that  the  pure  tone 
absorption  at  the  band  center  frequency  is  a  reasonable 
rough  approximation  to  the  absorption  encountered  by 
a  band  of  noise  analyzed  with  a  perfect  filter  provided 
the  band  is  narrow,  and  the  slopes  of  the  initial  source 
spectrum  and  the  air  absorption  curve  are  not  large. 

The  latter  constraint  implies  that  the  frequency  should 
be  low  and  the  propagation  path  short.  Specific  quanti¬ 
tative  guidelines  as  to  when  these  conditions  may  be 
met  will  be  given  later. 

B.  Varying  spectrum  slope 

If  the  spectral  density  P\{f)  of  the  source  has  a  vary¬ 
ing  slope  over  the  bandwidth  of  the  filter,  an  accurate 
computation  of  A,  for  this  case  requires  the  use  of 
some  type  of  spline  curve-fitting  function  for  estima- 
ting  at  each  of  the  Integration  frequencies  f,  over 
the  filter  bandwidth.  Several  methods  could  be  em¬ 
ployed  for  generating  such  spline  functions  including  a 
standard  cubic  spline  function  for  interpolating  between 
two  points.  However,  it  was  more  convenient  and  suf¬ 
ficiently  accurate  to  use  a  simple  log-linear  equation 
with  a  linearly  varying  elope  to  interpolate  between 
three  points— the  center  frequencies  of  three  adjacent 
filter  bands  centered  oo  any  one  band  of  interest.  Thus, 


referring  to  Fig.  2,  over  the  frequency  range  of  ±1  fil¬ 
ter  band  about  the  ith  band,  or  from to/,.,,  the 
spectrum  level  of  the  source  L,(f)  was  assumed  to  vary 

as 

!,(/)  -  £.,(/,)  + 10 [a,  +  6, (///,)]  logI0(///,) ,  dB  .  (21) 

where  L,(ft)  is  the  spectrum  level  at  the  center  fre¬ 
quency  (fc)  of  the  ith  band,  and  a, ,  b ,  are  the  constants 
determined  by  the  differences,  8,_  and  a,.,  between  the 
spectrum  levels  of  the  three  adjacent  bands  as  shown 
by  the  following  expressions: 

Bi  -  r’B,  J/10(r*-  l)log10r ,  (22a) 

*,«(«,♦- 6,  J/10(r»- l)log10r , 
and 

*<♦  *  ,)  —  £,(/,)  i 

6 

Note  that  when  8,.  “8,.,  then  6,  =0,  and  Eq.  (21)  would 
describe  the  spectrum  level  curve  of  a  linear  plot  of 
spectrum  level  versus  log  frequency  with  a  constant 
slope  of  10a,  dB/decade.  However,  the  type  of  smooth 
curve  depicted  in  Fig.  2  for  the  spline  interpolation 
method  is  readily  described  by  Eq.  (21)  for  nonzero 
values  of  8,.  This  more  general  approach  for  inter¬ 
polating  the  spectral  density  at  frequencies  between 
adjacent  band  center  frequencies  was  particularly  use¬ 
ful  when  dealing  with  practical  spectral  shapes  as  dis¬ 
cussed  later.  This  was  accomplished  by  using  the  vary¬ 
ing  slope  concept  specified  by  Eq.  (21)  to  define  the 
spectral  density  at  the  source,  P\,,  at  each >th  frequen¬ 
cy//  •• 

P'«-P*jj,tur'"",>"'' ,  (23) 

where  P^,  is  the  pressure  spectral  density  at  the  center 
frequency /,  of  the  ith  band,  and  «,,8,  are  the  slope  con¬ 
stants  derived  with  Eq.  (22)  from  the  estimated  spec¬ 
trum  levels  of  the  source  at  each  band  center  frequency. 
[The  same  techniques  outlined  here  to  compute  A,  for 
the  general  case  of  a  varying  slope  spectra  can  be  used 
to  determine  the  true  relationship  between  source  band 


dB, 

dB. 


(22b) 
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range  lor  the  filter  transmission. 

The  nonunifonn  transmission  response  of  the  filter 
can  be  based  on  the  actual  transmission  characteristics 
of  the  filter  set  for  which  the  filter  correction  factor 
A,  is  desired.  However,  a  conservative  and  more  gen¬ 
erally  applicable  approach  would  be  to  base  the  filter 
transmission  response  upon  standard  requirements  such 
as  specified  ANSI  SI.  11  — 1971. l*  This  approach  was 
adopted  for  this  paper.  An  empirical  expression  was 
employed  which  predicted  the  filter  transmission  re¬ 
sponse  curve  identified  in  Fig.  4  as  curve  (a),  the 
“practical”  class  HI  filter  assumed  for  this  analysis. 
Also  shown  in  Fig.  4  is  the  minimum  transmission  loss 
required  for  such  class  Cl  filters  by  the  ANSI  Standard 
(curve  b)  and  the  transmission  loss  of  a  typical  analog 
filter  (curve  c)  which  conforms  to  this  standard. 

The  "practical"  filter  response  closely  follows  that 
of  an  actual  filter  near  the  filter  band  edge  frequencies 
and  falls  slightly  below  the  ANSI  required  minimum 


FIG.  4.  Transmission  loss  curve  (a)  for  s  hypothetical  ''Prac¬ 
tical*'  j. -octave  band  (liter  compared  to  (b'  the  minimum  re- 
qui red  loss  according  to  ANSI  Standard  Si.  II  and  (c>  the  actual 
transmission  loss  of  a  typical  analog  filter. 


levels  (assuming  perfect  filters)  and  corresponding 
spectrum  levels  at  the  band  center  frequencies  without 
requiring  the  usual  assumption  of  a  fiat  spectrum  over 
the  filter  bandwidths.) 

Thus,  the  spectral  density  at  the  source  can  be  com¬ 
puted  at  any  frequency  (within  the  rar«e /,.,-/,♦,)  using 
Eq.  (23)  and  the  integration  of  Eq.  (1)  carried  out  using 
the  numerical  integration  techniques  specified  by  Eqs. 
(3)— (12). 

The  final  expression  for  A,  for  the  case  of  the  vary¬ 
ing  spectrum  slope  can  now  be  gj  ven  as 

ti'nl/J*  1  +  10  log  [£*,]- io  ,<*[£  flj  (24) 

where  Br/,B,j  are  the  powers  in  the  ;th  element  at  the 
receiver  and  source,  respectively,  as  defined  by  Eqs. 

(9)  and  (8).  The  corresponding  values  of  P ^  are  defined 
.by  Eq.  (23).  The  energy  sum  of  these  elemental  powers 
over  the  6  segments  defines  the  true  band  powers  of  the 
receiver  and  source,  respectively,  within  the  ideal  fil¬ 
ters. 

III.  DERIVATION  OF  A, 

When  a  real  filter  is  used,  the  same  numerical  inte¬ 
gration  techniques  Just  defined  for  the  computation  of 
A,  are  again  utilized.  However,  the  integra.  in  the  nu¬ 
merator  o I  the  second  term  of  Eq.  (1)  for  A,  is  modi¬ 
fied  In  Eq.  (2)  in  two  ways.  The  integration  limits  are 
extended  to  cover  the  practical  transmission  range  of 
the  actual  filter  and  the  nonuniform  power  transmission 
loss  of  the  filter  is  included  within  the  integral. 

Approximately  t  and  5  times  the  band  center  frequen¬ 
cy  define  reasonable  limits  for  the  practical  transmis¬ 
sion  range  of  a  filter.  However,  just  one  filter  band¬ 
width  below  or  above  the  filter  band  edge  frequencies 
can  often  be  used  as  effective  limits  for  the  frequency 


loss  at  frequencies  well  removed  from  the  filter  band 
edges.  This  “practical”  filter  response  was  generated 
by  slightly  modifying  the  expression  in  Ref.  16  for  the 
minimum  required  transmission  loss  to  obtain 

Tfi(/)  =  10Iogl0[A*a[C(///t)-(/c/Cf)n,  dB  (25) 

where,  for  J-octave  band  filters  (see  Appendix  E  of 
Ref.  6  for  corresponding  constants  for  octave  band  fil¬ 
ters) 

/t  =  8/13,  B  =  2547,  and  C=-10'l'*°,  for / * /,’ , 

A>1,B,  C  =  0,  for/;  </</;, 

A*  8  13,  B  >2547,  and  C  *  10J,'“  for  /*  /, 

and  and  10 l,Kfn  respectively,  the  band 

edge  frequencies  of  the  "practical”  filter  which  fall 
slightly  inside  the  nominal  band  edge  frequencies/,  and 
/,  of  an  ideal  i-octave  band  filter.  Note  that,  outside 
the  passband  of  this  practical  filter,  the  constant  B  in 
the  preceding  equation  differs  slightly  from  the  stand¬ 
ard  value  in  Ref.  16  in  order  that  the  transmission  loss 
will  be  zero  at  frequencies  /,*  and /,  and  also  equal  to 
4.5  dB  at  the  nominal  band  edge  frequencies/,  and/,, 
close  to  the  value  of  4.34  dB  recommended  by  Sepmeyer 
for  practical  filters.*’ 

The  final  expression  for  the  filter  correction  factor 
h r(R)  at  the  receiver  is  formed  by  expressing  the  con¬ 
tinuous  integrals  of  Eq.  (2)  in  the  form  of  summations 
ever  each  uf  they  discrete  elements  as  follows: 

A,(R)  =  101og^jB'r>  -lOlog^  Br/j,  dB.  (26) 

where  is  the  band  power  at  the  receiver  in  the/th 
integration  element  as  defined  by  Eq.  (9),  but  with  the 
pressure  spectral  density  P2  at  frequency  /,  replaced 
by  the  filtered  value  P^l O'"  <t>>  /l°.  The  summation  for 
the  first  term  above  is  carried  out  over  nb  segments— n 
bands  to  cover  the  effective  transmission  response 
range  of  the  fitter  and  A  segments  per  band.  The  sec¬ 
ond  summation  term  is  identical  to  (he  first  term  in 
Eq.  (24),  where  the  segment  powers  Brl  are  defined  by 
Eq.  (9)  and  are  summed  only  over  the  ft  segments  in  the 
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TIO.  5.  Correction  factor  A,  (SI  at  source 
for  various  Input  spectrum  alopea. 


pat  stand  of  the  Ideal  fitter. 

Equations  (10)  or  (11)  for  the  slope  parameter  m,  are 
also  modified  to  include  the  additional  change  in  the  in¬ 
put  spectrum  due  to  the  fitter  response.  For  example, 
£q.  (11)  becomes 

m,  •  U.l/yi>  -  TL(/Jtl)  -  L,(Jt) *  TL(/,))/(10/b) 

xlog10r .  (2D 

Typical  results  of  applying  this  process  are  illustrated 
in  the  next  section.  Note  that  if  B,t  is  used  instead  of 
B,t  In  Eq.  (26),  the  total  correction  A,  instead  of  only 
is  computed.  This  procedure  is  frequently  used 
when  the  spectral  density  at  the  source  la  known. 

IV.  RESULTS 

The  filter  error  at  the  aource  A,(S).  before  any  sub¬ 
stantial  additional  spectrum  shaping  has  occurred  due 


to  air  absorption,  is  normally  very  small  for  input 
spectra  with  a  constant  slope,  sb  illustrated  in  Fig.  5. 

In  fact,  for  {-octave  band  filters  which  equal  or  exceed 
ANSI  class  m  specifications,  the  filter  error  at  the 
source  A,(S)  can  probably  be  considered  completely 
negligible  for  most  source  spectra.  However,  at  a  dis¬ 
tant  receiver,  quite  a  different  picture  emerges.  Fig¬ 
ure  6  illustrates  the  value  of  the  filter  error  at  the  re¬ 
ceiver  A,(Jt)  for  typical  Jet  aircraft  spectra”  that  might 
be  observed  on  approach  or  takeoff  over  the  propaga¬ 
tion  path  lengths  indicated.  The  figure  is  drawn  to  il¬ 
lustrate  how  the  true  received  spectrum  in  such  cases 
should  actually  have  a  very  rapid  cutoff  at  frequencies 
above  2000-5000  kHz.  Thus,  when  A, (A)  is  not  account¬ 
ed  for  in  spectrum  analysis  of  distant  aircraft  sounds, 
high  frequency  levels  may  be  significantly  overstated. 

This  is  illustrated  more  directly  in  Fig.  7  which  il¬ 
lustrate!  the  absorption  attenuation  as  a  function  of 
propagation  path  length  for  the  same  typical  jet  aircraft 
takeoff  noise  spectra  used  in  Fig.  6.  One-third-octave 
band  levels  that  would  be  measured  (ignoring  spherical 
spreading  losses)  are  shown  for  four  frequencies  with 
ideal  or  practical  fitters  for  15°C  and  70%  relative 
humidity.  For  comparison,  the  excess  attenuation  at 
the  band  (geometric)  center  frequency  is  also  shown. 
Comparir^  the  latter  with  the  curves  for  the  ideal  filter 
illustrates  the  magnitude  of  the  ideal  fitter  correction 
A,.  Comparing  the  single  frequency  curve  with  the  top 
curve  in  each  group  illustrates  the  magnitude  of  the 
combined  correction  factor  Af  ♦  A r(R)  at  the  receiver. 
Clearly,  these  errors  become  quite  large  at  high  fre¬ 
quencies  and  large  distances.  However,  when  these 
errors  ore  large,  the  effected  band  levels  are  usually 
low  and  hence  weighted  overall  noise  levels  may  be  only 
slightly  influenced  by  the  errors  discussed  in  this 
paper. 


FKJ.  i.  Typical  variation  la  fitter  eorr action  factor  A,  Ut>  at 
receiver  for  eourco  epectra  from  aircraft  propagatod  over 

veriest  dMaaoea. 


Figure  6  illustrates  another  general  indication  of  the 
magnitude  of  the  combined  error  in  meaeured  band  lev- 
ele  A,  4  A,(Jt)  that  occur*  when  one  assume*  that  ab- 
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FIG.  7.  Exc ess  sttemuaMoo  due  to  atmospheric  absorption  of 
^-octave  baud  levels  for  typical  turbofan  aircraft  takeoff  spec¬ 
trum—  cotnpai ring  predicted  loss  at  band  center  frequency  and 
levels  that  would  be  meaaured  by  ideal  and  practical  filters. 

sorption  attenuation  of  a  band  level  equals  that  of  the 
band  center  frequency.  Similar  to  Fig.  3,  this  figure 
shows  contours  of  the  combined  errors  A  =  A;  +  &f(R)  as 
a  function  of  the  average  slope  SA  of  the  total  air  ab¬ 
sorption  and  assumes  the  source  spectrum  has  a  con¬ 
stant  slope  S,.  Ofl  the  lower  right  corner  of  Fig.  8,  the 


FIG.  #.  Curves  of  constant  band  correction  factor  (A  In  dB) 
In  terms  of  the  slops  of  the  source  spectrum  and  atmospheric 
attenuation  curve  for  an  ANSI  Class  HI  ^ -octave  band  filter  at 
288  K  and  ?0%  relative  humidity. 


FIG.  9.  Approximate  boundaries  of  frequency  and  propagation 
distance  for  which  total  band  attenuation  correction  exceeds 
0.S  dB.  An  approximate  -ule  of  thumb  is  given  by  the  dashed 
lines  which  indicate  that  the  total  band  correction  exceeds 
0. 5  dB  when  the  product  of  propagation  distance,  in  km.  and 
tbe  square  of  frequency,  in  kHz,  exceeds  4  or  10  for  octave  and 
X-octave  bands,  respectively. 


combined  (negative)  slope  of  the  air  absorption  and  the 
source  spectrum  can  equal  or  exceed  the  positive  slope 
of  tbe  lower  filter  side  band  and  the  numerical  integra¬ 
tion  does  not  converge.  However,  for  real  spectra,  the 
slope  of  the  spectra  will  eventually  decrease  sufficiently 
so  that  one  always  measures  a  finite  band  level. 

Figure  9  provides  an  approximate  guideline  as  to  the 
frequency-distance  regimes  where  the  combined  filter 
error  will  exceed  0.5  dB  for  the  type  of  "prac¬ 

tical”  filters  assumed  for  this  analysis. 

One  final  area  of  potential  application  of  these  correc¬ 
tion  factors  for  bands  of  noise  subject  to  significant  ex¬ 
cess  propagation  losses  is  concerned  with  normalizing 
measured  spectra  from  one  weather  condition  to  that  of 
another  weather  condition.  A  related  problem  is  cor 
cerned  with  the  evaluation  of  propagation  losses  for 
bands  of  noise  utilizing  a  layered  atmosphere  model. 

As  a  typical  example  of  the  former,  assume  that  a  spec¬ 
trum  Is  measured  when  the  atmosphere  is  homogeneous, 
with  a  temperature  of  25°C,  the  relative  humidity  equals 
70%  and  an  ANSI  class  E  J-octave  band  filter  is  used. 
For  reference  conditions  of  15*C,  70%  relative  humidity 
and  a  perfect  J-octave  band  filter,  at  4  kHz,  and  a  prop¬ 
agation  distance  of  400  m.  the  measured  band  would  be 
within  a  half  of  a  dB  of  that  which  would  have  been  mea¬ 
sured  under  reference  conditions.  We  can  conclude, 
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FIG.  10.  Application  of  correction 
(acton  for  computing  propagation  loan 
of  beads  of  nolaa. 


therefore,  that  the  effect  will  only  be  significant  when 
the  propagation  distance  is  very  large  (-km’s),  the  fre¬ 
quency  is  very  high  (-kHz),  or  the  actual  and  standard 
filters  vary  greatly.** 

Finally,  Fig.  10  summarizes  the  application  of  the 
two  correction  factors  discussed  in  this  paper.  As  out¬ 
lined  in  the  chart,  correction  for  different  propagation 
paths  (or  different  weather  conditions)  is  accounted  for 
primarily  by  changes  in  only  the  ideal  filter  correction 
factor  A,.  For  layered  atmosphere  models,  the  filter 
correction  factor  A,(A)  is  only  applied  once  at  the  re¬ 
ceiver.  This  correction  term  could,  in  fact,  change 
when  correcting  measured  data  at  a  receiver  from  one 
weather  condition  to  another  since  the  shape  of  the 
spectrum  entering  the  filter  might  be  different. 

The  method  outlined  in  this  paper  for  treating  the  in¬ 
fluence  of  the  spectrum  shape  on  atmospheric  absorp¬ 
tion  of  bands  of  noise  should  provide  reasonably  accu¬ 
rate  results  except  when  a  pure  tone  is  imposed  on  the 
spectrum.  In  the  absence  of  a  pure  tone,  the  slope  of 
the  spectrum  is  typically  easy  to  estimate  to  within  2  or 
3  dB  per  bandwidth  and  the  band  absorption  varies  only 
slowly  with  spectrum  slope.  When  a  pure  tone  is  pre¬ 
sent,  the  band  containing  that  tone  can  be  integrated  nu¬ 
merically  with  some  estimate  of  the  pure  tone  frequency 
within  the  band  while  the  remainder  of  the  noise  spec¬ 
trum  can  be  corrected  for  atmospheric  absorption  as 
outlined  here. 


V.  CONCLUSIONS 

The  sound  level  of  a  band  of  noise  examined  with  a  fil¬ 
ter  after  propagating  through  the  atmosphere  is  a  func¬ 
tion  of  the  atmospheric  conditions,  the  length  of  the 
propagation  path,  the  characteristics  of  the  filter,  and 
the  shape  of  the  initial  noise  spectrum  as  well  as  the 
source  intensity  or  pressure  level.  For  propagation  at 
high  frequencies  or  over  long  distances,  the  attenuation 
of  the  source  band  at  the  receiver  due  to  atmospheric 
absorption  is  not  accurately  given  by  the  attenuation  of 
the  band  center  frequency.  For  illustration,  it  is  con¬ 
venient  to  break  the  potential  error  into  two  parts,  al¬ 
though  in  actual  calculations,  it  is  desirable  to  compute 
the  combined  correction  The  first  part,  A/( 

defines  how  the  true  band  level  measured  with  an  ideal 
filter  differs  from  the  pure  tone  absorption  computed  at 
the  band  center  frequency.  The  other  factor,  A,,  ac¬ 
counts  for  the  finite  transmission  of  practical  filters 
outside  their  paasband.  For  long  distances  or  high  fte- 
quencies,  the  band  loss  can  be  much  smaller  (or  only 
slightly  larger)  than  the  pure  tone  absorption  loss.  Act¬ 
ual  filter  characteristics  enhance  this  effect  making  it 
difficult  if  not  impossible,  in  extreme  cases,  to  deduce 
the  true  source  spectrum  from  values  observed  at  a 
distance  receiver.  Correcting  a  receiver  spectrum  to 
standard  conditions  is  more  complex  but  typically  may 
involve  only  a  small  correction  for  band  filter  effects. 

In  view  ot  the  rather  large  corrections  to  be  applied 
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to  measured  band  levels  to  account  (or  atmospheric  ab¬ 
sorption  in  extreme  cases,  estimates  of  this  quantity 
should  be  made  while  planning  a  series  of  measurements. 
The  propagation  correction  can  be  reduced  by  reducing 
propagation  distance  and  selecting  narrow -band  filters. 
The  fact  that  the  band  correction  is  largest  (and  some¬ 
times  not  well  defined)  for  high  frequencies  casts  doubt 
on  spectral  measurements  made  of  distant  signals  such 
as  aircraft  noise  or  shock  waves  with  significant  high 
frequency  components  at  their  source. 

This  paper  has  dealt  with  the  correction  of  a  source 
spectrum.  In  the  real  world,  one  has  only  a  received 
spectrum.  To  follow  the  reverse  path  will  require  an 
iterative  procedure  until  a  source  spectrum  is  found 
such  that  when  the  spectrum  is  corrected  for  atmo¬ 
spheric  and  filter  effects,  the  measured  received  spec¬ 
trum  results.  This  procedure  is  now  being  investigated 
with  particular  attention  to  the  question  of  whether  or 
not  the  deduced  source  spectrum  is  unique.  These  re¬ 
sults  will  be  reported  in  a  future  paper. 
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APPENDIX  C 


ANALYSIS  OF  POWER  TRANSMISSION  THROUGH  A 
PRACTICAL  FRACTIONAL-OCTAVE  BAND  FILTER 

1.0  INTRODUCTION 

The  filter  error  is  defined,  in  the  main  body  of  the  text,  os  the  difference 

in  decibels  between  the  power  actually  transmitted  by  a  practical  filter  and  an  ideal 

filter.  This  problem  is  treated  in  detail  in  this  appendix  for  the  case  of  various  constant 

slope  input  signals.  The  analysis  involves  numerical  integration  over  the  effective 

transmission  range  of  a  mathematical  approximation  to  the  real  transmission  response  of 

an  actual  filter.  This  approximation  closely  matches  the  actual  transmission  response  of 

current  technology  6-pole  Butterworth  filters  for  frequencies  extending  roughly  between 

the  20  dB  down  points  of  the  filter  response  and  then  tends  to  provide  a  conservative 

20 

approximation  to  the  filter  response,  within  ANSI  standard  specifications,  for  lower 
portions  of  the  filter  skirts.  It  is  possible,  with  a  digital  filter  technology,  to  provide 
much  sharper  filter  skirts  than  is  provided  by  a  6-pole  Butterworth  filter.  (The  latter 
has  a  roll-off  rate  approaching  about  70  to  80  dB/octave  in  the  steepest  portions  of  the 
filter  response.)  However,  current  commercial  practice  is  often  limited  to  simulating, 
digitally,  no  better  than  6-pole  Butterworth  filter  response  characteristics. 

It  should  be  pointed  out  that  the  specification  for  the  minimum  transmission 
response  characteristics  of  one-third  octave  band  filters  to  be  used  for  aircraft  noise 
certification  are  specified  by  the  International  Standard  IEC  Publication  225 •  This 
is  less  stringent  in  roll  off  characteristics  than  that  specified  in  the  comparable 
American  Standard  SI.  11  for  frequencies  where  the  attenuation  is  -20  dB  or 
more  below  the  pass  band  response.  (These  response  characteristics  ore  compared  in 
Figure  1  of  this  appendix.)  Thus,  the  results  outlined  in  this  appendix  would  be 
somewhat  optimistic  for  spectrum  analysis  filter  sets  which  met  the  IEC  Standard 
but  not  the  ANSI  Standard.  In  fact,  however,  most  commercially  available  filter 
systems  meet  both  standards,  so  that  these  results  are  actually  slightly  pessimistic  for 
spectra  with  large  negative  slopes. 
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2.0  INTEGRATED  POWER  TRANSMISSION 


log  T(f) 

loTwf) 


2.1 


where 


and 


where 


Over  any  bandwidth  -  f^,  the  total 
power  B  transmitted  by  a  filter  with 
a  power  transmissibility  function  T(f) 
driven  by  a  broadband  noise  with  a 
spectral  density  W(f),  as  indicated  in 
the  sketch,  is  given  by 


B  *  /  T(f)  W(f)  •  df  (1) 


Constant  Slope  Cose  —  Ideal  Filter 

Let  the  input  spectrum  over  any  band  fj  to  f^  be  given  in  the  form 

W(f)  =  W1  *(f/f,)m  (2) 


W]  =  W(f)  at  f  =  f 

the  filter  response  over  the  same  band  be  given  by 

T(f)=T1*(f/f1)n  (3) 


T,  «T(f)  at  f  =  fj 
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Then,  for  this  constant  slope  case,  the  power  (B)  transmitted  through  the  filter  in  this 
band  can  be  defined  in  closed  form  by  the  integral 


V,  fr— 1 

^  m+n  <1  ^  m+n  L  m  +  n  +  J  J 


or  B  =  W17l  f1  [W"  +  "  +  1  '  >] /(m  *  "  +  >>  (*) 

2.2  General  Cose 

For  the  more  general  case,  such  as  for  the  power  transmitted  through  an  actual 
filter,  the  filter  response  function  is  definable  in  closed  form  as  a  smooth  curve  but  the 
result  is  not  readily  integrable.  A  piece-wise  integration  is  therefore  necessary  using 
the  above  results  for  each  segment. 

2.2.1  ANSI  Filter  Response 

The  power  transmission  curve  for  ANSI  filters,  considering  the  minimum  trans¬ 
mission  loss  boundary,  can  be  defined  by  an  equation  of  the  form  (see  Appendix  A) 


TL  *  10  logI0CT(f)]  =  -l0log10 


,  dB 


(5) 


where  a,  b,  k  are  constants  for  each  portion  of  the  filter  response  curve  and  f.  = 
geometrical  mean  frequency  of  the  filter  band. 

There  are  three  regions  for  the  filter  response  curves  as  indicated  in  the  sketch 

below. 
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P  the  "zero  loss"  pass  band,  F.j  lo  f. ^  where  F.^/f.  =  f./f. ^ 

S  the  "stop"  bands,  Pj2  to  F^  and  to  f.j  where  *  f./fj2 

L  the  limiting  stop  bands,  <  P.2  and  > 


General  Shape  of  Frequency  Response  of  Filter  with 
Minimum  Transmission  Loss  Allowed  by  ANSI  Specification 


The  following  table  summarizes  the  constants  a,  b,  k,  and  the  frequencies,  relative 
to  f.,  for  ANSI  Class  III  one-third  octave  band  and  Class  II  octave  band  filten  for  each  of 

i 

these  three  regions. 

Table  1 


Constants  for  Minimum  Transmission  Loss  Curve  for  ANSI  Filters* 


*See  Appendix  C-1. 
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Figures  1  and  2  show  the  frequency  response  curves  corresponding  to  these  minimum 
attenuation  requirements.  Response  curves  for  typical  actual  filters  are  also  illustrated. 

For  conservatism,  it  is  desirable  to  use  a  response  curve  just  slightly  better  than  the  mini¬ 
mum  attenuation  specification.  This  provides  an  approximation  to  a  real  filter  without 
the  use  of  any  one  existing  model  and  also  insures  a  conservative  result  for  the  filter 
effect. 

2.2.2  "Practical"  Filter  Response 

To  approximate  the  "practical  filter, ”  it  was  convenient  to  simply  modify  Equa¬ 
tion  (5)  by  multiplying  the  frequency  f  by  a  scaling  factor  c  which  has  the  following 
values.  (The  exact  values,  used  for  c,  equal  to  powers  of  10,  are  defined  along  with  the 
approximately  corresponding  values  as  powers  of  2.) 


)0l/60  „  2t/18 

r  f>f. 

For  one-third  octave  band  filters,  e  =| 

10-|/6°.  2-'/l8 

,  f<  f. 

,o'/20-2,/6 

,f>f. 

For  octave  band  filters,  c  - 

l10-l/20-2',/6 

,  f<  f. 

The  "zero  transmission  loss"  frequency  bandwidths  decrease  accordingly  so  that 

the  frequencies  F. .  and  f..  become 
•  I  il 

For  l/3rd  octave  band  filters,  P  =  2^®  •  2  ^  *  f.  *  2  ’  f  *■  10  f. , 

il  iii 

and  f-10,/30f.  -  1.0798  f. 

il  i  i 

For  octave  band  filters,  F.,  -  21/6  *  2“1/2  •  f.  »  21/3  f.  -  10_1/10  •  f., 
and  f  =10,/10f.  -  1.2589  f. 

il  i  i 
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For  the  "practical"  filter,  it  was  assumed  that  the  filter  response  in  the  limiting 
stop  bands  (region  L)  could  be  extrapolated  from  the  response  in  the  stop  bands  (region 
S)  by  using  the  constants  in  Equation  (5)  appropriate  for  this  region  to  compute  the 
response  for  the  limiting  stop  bands  (region  L).  This  extrapolation  was  carried  out  to 
upper  and  lower  limiting  frequencies  which  were  set  approximately  equal  to 
10  f.  and  fj/10.  This  arbitrary  procedure  was  mode  to  avoid  unrealistic  power  transmission 
of  the  filter  in  the  limiting  stop  bands.  In  actual  practice,  the  noise  floor  of  a  spectrum 
analysis  system  will  often  be  comparable  to  the  power  transmitted  by  a  filter  at  these 
extreme  frequency  regions  on  either  side  of  the  pass  band  of  a  filter. 

In  summary,  then,  the  "practical"  filter  response  is  computed  with  the  following 
equation  and  parameters  indicated  in  Table  2  for  each  filter  type  and  frequency  range. 


Table  2 

Constants  for  Computing  Transmission  Loss  of  "Practical"  Filters  Which 
Approximate  the  Minimum  Transmission  Loss  of  ANSI  Filters 


•The  values  of  b  differ  slightly  from  those  in  Table  1  to  achieve  on  attenuation  of  0  dB  for  these  practical 
Filters  at  the  limiting  frequencies  f'.j  and  f,j  • 


The  resulting  transmission  loss  curves  ore  also  shown  in  Figures  1  and  2  compared 
to  the  response  curves  for  typical  filters  and  the  ANSI  specifications  for  minimum  loss. 

The  same  information  is  also  shown  in  Figures  3  and  4  on  linear  scales  for  the  transmission 
coefficient  T(f)=  log^  ' (T L/l 0)  and  frequency.  On  this  type  of  plot,  it  is  easy  to 
recognize  the  relative  power  transmitted  for  each  type -df  filter  according  to  the  areas 
under  the  corresponding  response  curve. 

2.2.3  Integration  of  Filtered  Spectrum 

To  determine  the  power  transmitted  by  this  "practical"  filter  for  an  arbitrary 
input  spectrum  W(f),  the  entire  frequency  range  will  be  broken  up  into  narrow  constant- 
percentage  frequency  bands  with  a  constant  ratio  between  the  upper  and  lower  frequencies 
for  each  segment.  This  frequency  ratio  r  will  be  equal  to 

J  A> 


r=  10' 


(7) 


where  b  =  an  even  integer. 


Within  plus  or  minus  one  filter  bandwidth  on  each  side  of  the  nominal  band- 
edge  frequencies  of  a  filter,  the  constant  b  is  set  equal  to  a  value  of  240  for  one-third 
octave  band  filters  and  60  for  octave  band  filters.*.  Outside  this  critical  range,  wider 
integration  segments  are  satisfactory  so  that  b  is  decreased  to  a  minimum  value  of  20. 
Thus,  for  an  octave  band  filter,  power  transmission  was  computed  over  J/60th  decade 
(~  1/1 8th  octave)  segments  for  frequencies  between  (1/(2  Vi”))f.  and  ( 2  *  V  2  )  •  f., 
(i.e.,  db  1  octave  on  other  side  of  the  pass  band)  and  in  l/20th  decade  (~  l/6th  octave) 
segments  outside  this  frequency  range.  It  was  convenient  to  use  integration  segments 
with  band-edge  frequencies  defined  in  terms  of  powers  of  10  instead  of  powers  of  2  in 
order  to  provide  a  better  match  to  one-third  octave  band  filter  center  frequencies. 


4 An  initial  evaluation  of  integration  of  an  arbitrary  spectrum  over  an  octave  or  one- 
third  octove  bandwidth  indicated  that  numerical  integration  with  segments  narrower 
then  1 /240th  of  o  dec  ode  (or  l/72nd  of  an  octave)  changed  the  resulting  integrated 
power  by  lees  than  0. 01  dB . 
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Figure  3.  Relative  Power  Transmission  for  One-Third  Octave  Band  Filter 
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f/fj,  Frequency/Bond  C«n**r  Frequency 


Figure  4.  Relative  Power  Trorsmiision  far  Octave  Band  Filter* 


As  illustrated  in  Figure  5,  it  was  assumed  that,  over  each  segment,  the  input 
spectrum  and  filter  transmission  response  could  be  described  by  simple  power  law 
expressions,  as  given  by  Equations  (2)  and  (3),  corresponding  to  constant  slopes  of  the 
spectral  density  W(f)  and  the  transmission  function  T(f)  versus  frequency  when  these  are 
plotted  on  logarithmic  scales. 


Figure  5.  Illustration  of  Integration  of  Filtered  Spectrum  W'(f)  over  Frequencies 
f.  to  f.  +  j  Using  Straight  Line  Segments  for  Input  Spectrum  Shape  W(f) 

and  Filter  Response  T(f) 


For  the  j  segment,  the  power  transmitted  (B.)  can  be  written  down,  by 
inspection  of  Equation  (4)  as 
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8j  -  W(f.)  T(f.)  f.  [(f.  +  ,/fjf  *  "  +  '  -  »]/(•>  *  "  +  » 


where 

W(f.)  =  Input  spectral  density  at  frequency  f. 

T(f.)  =  Transmissibility  of  filter  at  frequency  f, 

m  =  Slope  of  input  spectral  density  from  f .  to  f  +  j 
n  *  Slope  of  filter  response  from  f.  to  f^  +  j 

For  a  constant  slope  S  of  the  input  bond  level  spectrum  in  dB  per  bandwidth, 
it  can  be  shown  that  the  slope  exponent  (m)  for  the  input  spectral  density  is 

m  «  (S/10  log  r')-l  (9) 

1/3 

where  r'  =  ratio  of  center  frequencies  of  adjacent  bands  (2  or  2  '  for 
octave  and  one-third  octave  filters/  respectively) 

The  slope  exponent  (n)  for  the  transmissibility  function  T(f)  can  be  given  by 

,°«io(V*i/,i) 


where 


.nft*i)-aft) 

NWi) 


Tiff.  =  the  transmission  loss  of  the  practical  filter  as  given  by 
'  +  '  Equation  (6)  and  Table  2  for  f  =  f  ^ 


•For  the  special  case  of  m  +  n  +  l-~-0,  it  con  be  shown  that  Equation  (8)  reduces 

to  B.—»  W(f.)T(f.)  f.  In  [f.../f.] 

J  I  I  J  •  J  +  1  J 
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In  the  case  of  a  more  general  input  spectrum  which  does  not  have  a  constant 
slope,  then  Equation  (10)  can  be  applied  to  obtain  the  sum  of  the  two  exponents  (n  +  m) 
by  changing  the  numerator  to  include  both  the  input  spectrum  W(f)  and  filter  response 
T(f)  as  follows 


n  +  m  = 


So[w<fi 


+  1)-T(f.  +  |)/W(f.)-T(f.)] 

HoW? 


(H) 


It  will  be  convenient,  in  all  cases,  to  present  the  resulting  final  integrated 
power  through  the  filter  in  a  normalized  form  by  dividing  the  true  power  output  by  a 
reference  power  output  for  an  ideal  filter  with  a  constant  input  spectral  density  equal  to 
the  value  W.  at  the  geometric  mean  frequency  (f.)  of  the  band.  Thus,  summing  over  all 
j  segments,  this  relative  power  transmission  (Br)  will  be 


B 

r 


X*: 


Wi  fi 


(12) 


and  0  = 


/  2  ^  for  octave  band  filters 


(2'^  -  1)  /  2}^  for  one-third  octave  band  filters 


|  The  bandwidth  of  the  ideal  filter  relative  to  its  geometric  mean 
center  frequency  f. 


B.  is  the  power  transmitted  in  the  segment  as  defined  by  Equation  (8)  using 
Equations  (9),  (10),  or  (1 1)  to  define  the  exponents  n  and  m,  and  Equation(6)  and  Table  2 
to  define  the  filter  response  TL  -  10  log^  [T(f)]  .  The  input  spectrum  W(f)  is  specified 
for  the  case  of  a  constant  spectrum  slope  by 


w(f)  =  w.  .•  (f/f.  )m 


(13) 


For  this  case,  combining  Equations  (8),  (12),  and  (13),  the  relative  power  transmission  is 
the  summation 


•r-Z  T('i)  •  (</<;T  +  '  [(fi  ♦  +  "  *  '  -  ']/["  +  "  + ']  e  <“> 

and  the  summation  is  carriedout  over  all  the  j  segments*  necessary  to  encompass  the  entire 
frequency  range  for  f.  from  f./lO  to  10  *  f. . 

This  filter  response  analysis  has  been  carried  out  for  the  approximations  to  ANSI 
octave  and  one-third  octave  filters  discussed  earlier  and  for  a  range  of  input  spectra  with 
a  constant  slope  of  the  true  band  level  varying  from  -24  dB  to  +12  dB  per  octave.  (The 
corresponding  slope  of  the  spectrum  level  is  just  3  dB  per  octave  less.)  The  results  are 
presented  in  Table  3. 

In  order  to  provide  a  better  picture  of  the  power  transmission  through  the 
"practical  "  filters,  the  summation  involved  in  Equation  (14)  has  been  carried  out  over 
discrete  frequency  intervals  corresponding  to  nominal  bandwidths  of  the  filter  itself.  Thus, 
in  Table  3a,  reading  down  the  second  column,  the  power  transmission  is  shown  in  each  of 
7  octave  bands/  numbered  -3,  -2,  -1,0,  1,2,  and  3  including  the  center  pass  band  (band  0), 
(covering  a  range  of  *1  decade  about  the  nominal  center  frequency  of  the  filter)  for  an 
input  spectrum  with  a  band  level  slope  of  -24  dB  per  octave.  The  bottom  row  of  the  table 
gives  the  total  sum  from  all  seven  of  the  band  segments  which  corresponds  to  the  total  power 
output  that  would  be  measured  by  this  'practical  "  filter.  As  discussed  earlier,  the  numbers 
in  these  tables  are  actually  the  power  transmitted  relative  to  the  power  transmitted  by  an 
ideal  filter  with  a  white  noise  input  having  the  same  spectral  density  as  the  actual  spectra 
density  at  the  band  center  frequency  (fj).  The  bandwidth  of  this  ideal  filter  is,  of  course, 
the  nominal  bandwidth  of  the  actual  filter  (i.e.,  fj/^2  for  on  octave  filter). 

A  brief  examination  of  the  tables  show  that, as  one  can  expect,  they  exhibit  odd 
symmetry  about  the  intersection  of  the  0  pass  band  row,  and  0  db/octave  column.  The 
relative  power  actual  transmitted  through  the  pass  band  is  about  0.6  dB  less  than  for  the 
ideal  filter  due  to  the  rounding  of  the  response  curve  at  the  band  edge  of  the  practical 
filter.  However,  for  a  white  noise  input  (band  slope  equal  +3  dB/octave),  the  total 

*See  first  paragraph  in  Section  2.2.3,  page  7,  for  definition  of  number  of  segments 
used  for  integration. 
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Toble  3 

A.  Power  Transmission  for  ANSI  Closs  II  Octave  Bond  Filter  -  dB 


REFERENCE  IS  0  SB  THRU  IDEAL  PASS  BOMB  FILTER  (BAND  0>  FOR  WHITE  NO ICE  INPUT 
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B.  Power  Transmission  for  ANSI  Class  III  1/3  Octave  Band  Filter  -  dB 
REFERENCE  IS  0  DB  THRU  IDEAL  PASS  BAND  FILTER  (BAND  0>  FOR  UNITE  NOISE  INrUT 
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power  transmitted  for  either  the  octave  or  one-third  octave  practical  filter  is  only  about 
0.05  to  0.06  dB  below  that  of  the  ideal  filter. 

The  principal  result  obtained  from  the  analysis  is  to  clearly  show  the  substantial 
portions  of  the  total  power  output  from  a  "practical  "  filter  for  on  input  with  a  constant 
spectrum  slope  that  is  transmitted  by  the  filter  skirts,  particularly  the  portions  of  the  filter 
response  within  ±1  filter  band  of  the  pass  band. 

Based,  in  part,  upon  these  results,  the  spectrum  iteration  method  defined  in 
Section  3  of  this  report,  was  developed  to  moke  it  possible,  during  data  analysis,  to 
account  for  filter  response  characteristics  directlv  so  that  it  would  not  be  necessary  to 
apply,  operationally,  the  filter  error  term  Aj.  which  varies  directly,  in  a  complex  way, 
with  the  signal  spectrum  shape,  and  hence,  indirectly,  with  propagation  path  and  weather 
conditions. 

One  other  point  should  be  clear  from  this  appendix  -  the  exact  results  obtained 
are  uniquely  defined  by  the  filter  response  characteristics  assumed.  As  shown  earlier,  the 
approximation  used  in  this  report  is  considered  a  close  approximation.  However,  before 
any  routine  application  of  the  spectrum  iteration  technique  could  be  applied  in  industry, 
further  examination  of  this  filter  response  description  might  be  desirdble. 
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APPENDIX  C-l 


MINIMUM  AND  MAXIMUM  TRANSMISSION  LOSS  (TL) 
FOR  ANSI  ONE-THIRD  AND  FULL  OCTAVE  BAND  FILTERS 


1.  ANSI  Third-Octave  Band  Filter?  -  Class  III* 

(1)  For  <  j-  <  ~  (Pass  Band),  die  transmission  loss  (TL)  is  specified  by: 


Oi  TL  <  10  log 


5 

7 


,  dB 


(2)  For  ^  ^  ^r-  S  2  ^  or  2  ^  s  y—  £  5  (Stop  Bands) 

TL  >  10  log  £j|-  +  250o(j-  -  ^  J  ,  dB  -  0@  £-  *  2  1/6or  2_1/6 

(3)  For  y-  <  j  or  y-  >  5,  (Limiting  Stop  Bands) 

i  i 

TL  >  75  dB 


2.  ANSI  Octave  Band  Filters  -  Class  II* 

3  f  4 

(1)  For^*^-*j  ~  Pass  Band  Region 


0.n«i0h«[|  *t(^-t)] 


,  dB 


‘Americas  National  Standard  Specification  for  Octave,  Half-Octave,  and  Third- 
Octave  Band  Filter  Sets.  ANSI  SI .  1 1- 1966(R-1971 ). 
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